Technical  Report  1456 

October  1991 

User’s  Guide  for  the 
VTRPE  (Variable  Terrain 
Radio  Parabolic  Equation) 
Computer  Model 


Frank  J.  Ryan 


Approver!  for  public  release;  distribution  is  unlimited. 


NAVAL  OCEAN  SYSTEMS  CENTER 

San  Diego,  California  92152-5000 


J.  D.  FONTANA,  CAPT,  USN  R.  T.  SHEARER 

Commander  Technical  Director 


ADMINISTRATIVE  INFORMATION 

The  computer  code  described  in  this  report  was  developed  through  funding 
provided  by  the  Joint  Electronic  Warfare  Center  (JEWC)  under  MIPR  90-0069,  and 
by  the  Office  of  Naval  Technology  (ONT)  under  the  ONT  block  program  N02C  at 
the  Naval  Ocean  Systems  Center  (NOSC). 

Released  by  Under  authority  of 

M.  F.  Morrison,  Head  J.  H.  Richter,  Head 

Acoustic  Branch  Ocean  and  Atmospheric 

Sciences  Division 


EXECUTIVE  SUMMARY 


This  report  is  a  user's  guide  to  the  VTRPE  (variable  terrain  radio  parabolic  equation) 
computer  model.  It  is  designed  to  provide  the  reader  with  a  summary  of  the  physics 
and  numerical  methods  used  in  the  VTRPE  model,  along  with  detailed  instructions  on  the 
model’s  use  and  operation. 

The  VTRPE  computer  program  is  a  range-dependent,  tropospheric  microwave  propaga¬ 
tion  model  that  is  based  upon  the  split-step  Fourier  parabolic  wave  equation  algorithm. 
The  nominal  applicable  frequency  range  of  the  model  is  YHF  to  K-band.  The  VTRPE 
program  is  able  to  make  predictions  for  microwave  propagation  over  both  land  and  wa¬ 
ter.  The  VTRPE  code  is  a  full-wave  propagation  model  that  solves  the  electromagnetic 
wave  equations  for  the  complex  electric  and  magnetic  radiation  fields.  The  model  accounts 
for  the  effects  of  nonuniform  atmospheric  refractivity  fields,  variable  surface  terrain,  and 
varying  surface  dielectric  properties  on  microwave  propagation. 

Multiple  refractivity  vs.  altitude  profiles  may  be  input,  and  the  VTRPE  code  will  inter¬ 
nally  compute  a  smooth,  continuous  rrefractivity  field  in  range  and  altitude.  The  surface 
dielectric  properties  may  ’ -e  input  directly  into  the  model  or  computed  internally  using 
semi-empirical  algorithm'.  For  overland  propagation,  the  terrain  elevation  profile  is  user 
selectable. 

A  wide  variety  of  transmitter  characteristics  are  encompassed,  including 

•  horizontal  and  vertical  polarization. 

•  variable  mainlobe  vertical  beamwidth  and  tilt. 

•  arbitrary  transmitter/receiver  geometries,  and 

•  analytic  antenna  radiation  patterns. 

The  highly  efficient  split-step  Fourier  parabolic  equation  algorithm  is  used  to  calculate 
the  antenna  pattern  propagation  factor,  which  is  a  critical  component  of  the  radar  trans¬ 
mission  equation,  that  is  used  in  the  performance  prediction  and  analysis  of  radar  and 
communication  systems.  Optimized  fast  Fourier  transform  (FFT)  algorithms  are  used  in 
the  code  for  high  computational  efficiency.  The  user  has  control  over  both  the  FFT  size 
and  horizontal  range  step  size  used  in  the  calculation.  Many  model  inputs  have  internal 
defaults,  thus  minimizing  the  amount  of  input  data  required  for  execution. 

The  code  is  higlily  portable,  being  written  in  AN'SI-77  FORTRAN  with  MILSPEC- 
1753  FORTRAN  language  extensions.  Modern  software  structures  and  modular  program 
design  are  used  throughout.  The  VTRPE  program  is  currently  configured  to  run  under 
the  UNIX  operating  system  on  SUN  minicomputers  and  CONVEX  supercomputers,  and 
under  MS-DOS  on  8038G/8048G- based  PC's. 
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§1.0  Introduction 


This  report  is  designed  to  be  a  guide  to  the  use  and  operation  of  the  VTRPE  (vari¬ 
able  terrain  radio  parabolic  equation)  computer  code.  The  VTRPE  code  is  a  microwave 
propagation  model  designed  for  use  in  range- dependent  environments.  It  provides  a  full- 
wave  solution  (i.e„  the  electromagnetic  wave  equations  are  solved  for  the  complex  fields) 
to  the  problem  of  microwave  propagation  in  a  range-dependent  environment  having  vari¬ 
able  surface  boundary  conditions.  The  model  is  based  upon  the  parabolic  equation  (PE) 
or  paraxial  approximation  to  the  electromagnetic  wave  equation,  and  uses  the  split-step 
Fourier  PE  algorithm.  The  VTRPE  model  is  an  extension  of  the  smooth  surface  RPE  (radio 
parabolic  equation)  model  developed  by  the  author  to  predict  microwave  propagation  in 
marine  environments.1  The  major  difference  between  the  VTRPE  and  RPE  codes  is  in  the 
ability  of  the  former  to  handle  complicated  finite  conductivity  boundary  conditions  and 
irregular  surface  terrain. 

The  primary  use  of  the  VTRPE  program  is  to  predict  tropospheric  microwave  propa¬ 
gation  in  range-dependent  environments.  This  information  is  used  in  the  analysis  and 
performance  prediction  of  radar  and  communication  systems  operating  in  the  microwave 
region  of  the  spectrum.  The  frequency  range  of  interest  is  YHF  to  Ivu-band  (~  50  MHz  - 
IS  GHz).  At  these  frequencies,  spatial  variations  in  the  tropospheric  refract ivity  field  and 
the  earth's  surface  dielectric  properties  can  have  major  effects  on  radar  signal  propagation. 
In  fact,  substantial  differences  may  exist  between  radar  propagation  in  range- dependent 
environments  and  the  classic  “standard  atmosphere”  that  is  typically  assumed  by  radar 
system  designers. 

In  practice,  the  actual  performance  of  a  radar  in  the  atmosphere  is  often  quite  different 
from  that  predicted  based  upon  free-space  propagation  conditions.  Free-space  ranges  are 
often  several  orders  of  magnitude  different  from  observed  detection  ranges  in  the  atmo¬ 
sphere.  The  reason  for  this  discrepancy  is  threefold: 

1.  The  earth's  surface  is  a  finite  conductor  and  scatters  (reflects)  incident  energy  in 
various  directions,  leading  to  complicated  spatial  interference  patterns. 

2.  The  curved  earth  casts  a  shadow  giving  rise  to  diffraction  phenomena. 

3.  Inhomogeneities  in  the  atmospheric  index  of  refraction  cause  significant  refraction 
or  bending  of  radio  wave  energy. 

The  VTRPE  model  incorporates  full-wave  propagation  physics  to  properly  account  for  these 
"anomalous "  propagation  effects. 

1.1  Antenna  Patterns 

In  analyzing  these  various  phenomena,  it  proves  useful  to  separate  those  system  param¬ 
eters  not  influenced  by  the  environment,  for  example  the  antenna  radiating  characteristics, 
from  propagation  effects  that  are  environmentally  influenced,  such  as  ducting. 

The  radiating  characteristics  of  an  antenna  are  specified  in  terms  of  the  antenna  radia¬ 
tion  pattern  function  f{6.  o).  where  (6.  o)  are  the  zenith  and  azimuthal  angles,  respectively, 
of  a  spherical  coordinate  system  centered  at  the  antenna,  with  polar  axis  pointed  in  the 
direction  of  maximum  transmission.  The  antenna  radiation  pattern  function  is  defined  to 
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be  the  ratio  of  electric  (or  magnetic)  field  strength.  E(0.  0).  radiated  in  the  direction  (0, 0) 
to  the  peak  transmitted  field  strength,  E0  : 

a) 

The  antenna  pattern  function  /  is  related  to  the  time-averaged  Poynting  vector  5  of 
the  electromagnetic  wave  field  by 

S(8,<i>)  =  |/(0,  c>)|2S0. 

where  S0  is  the  energy  flow  per  unit  area  corresponding  to  the  peak  field  E().  In  general, 
the  antenna  radiation  pattern  function  /  is  a  complex  valued  quantity. 

Closely  related  to  the  antenna  pattern  function  is  the  antenna  gain.  The  antenna  gain 
G  is  expressed  in  terms  of  the  antenna  radiation  pattern  function  /  as 

r  4t 

I/<#.o)|2<B1' 

Note  that  the  antenna  radiation  pattern  function  /  and  the  corresponding  antenna  gain  G 
are  conventionally  defined  with  respect  to  free-space  propagation  conditions,  and  therefore 
do  not  include  anv  environmental  effects. 


1.2  Radar  Transmission  Equation 


To  systematically  incorporate  propagation  effects  and  antenna  characteristics  in  radar 
system  performance  calculations.  the  radar  transmission  equation  is  oiten  employed.  Fol¬ 
lowing  Kerr.2  define  the  one-way  generalized  transmission  equation,  which  relates  the 
power  received  by  an  omnidirectional  receiver  at  a  point  in  space.  Pr ( r ) .  to  the  power 
emitted  from  a  transmitting  antenna.  Pt.  by 

PM)  ..  r  f  Ft  l2  ,0, 

~pr-G'[ 2M?J  •  {2) 

and  the  corresponding  two-way  generalized  transmission  equation  for  a  monostatic  radar 
(i.e..  the  transmitter  and  receiver  are  colocated): 


where 


=  G'G^[ 


Gj  =  transmitting  antenna  power  gam. 
Gr  =  receiving  antenna  power  gain. 

27T 

k'n  =  vacuum  wavenumber  =  . 


a  =  target  radar  cross  section. 

R  =  radar-to-target  range. 

Fr  =  pattern  propagation  factor  for  target-to-receiver  path. 

Ff  =  pattern  propagation  factor  for  transmitter-to-target  path. 
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Environmental  effects  are  included  in  the  generalized  transmission  equation  via  the 
pattern  propagation  factor  F.  The  pattern  propagation  factor  accounts  for  the  fact  that 
(1)  the  receiver  may  not  be  in  the  mainlobe  beam  maximum  of  the  vertical-plane  antenna 
pattern,  and  (2)  non- free- space  wave  propagation  may  occur. 

The  pattern  propagation  factor  F  is  defined  as  the  ratio  of  the  field  magnitude  at  a 
point  in  space,  E(r).  to  the  magnitude  of  the  field  at  the  same  point  but  under  free-space 
propagation  conditions.  E0(r): 

F  =  .  (4) 

Eo(r) 

This  definition  of  the  pattern  propagation  factor  assumes  that  the  transmitting  antenna 
is  aligned  with  its  maximum  response  axis  pointed  directly  at  the  observation  point.  The 
pattern  propagation  factor  F  is  the  component  of  the  generalized  transmission  equation. 
Eq.  (2)  or  Eq.  (3).  that  includes  antenna  directivity  and  medium  propagation  effects.  It  is 
the  fundamental  quantity  to  be  computed  by  the  VTRPE  model.  The  other  parameters  in 
the  radar  transmission  equation  are  independent  of  the  environment. 

Because  of  the  large  dynamic  ranges  of  the  various  terms  in  the  radar  transmission 
equation,  it  is  customary  to  work  in  dB-space.  Thus,  define  the  dimensionless  variables 
PF  and  PL  by 

PF  =  20  log  |F|,  (5) 

PL  =  201og(2 k0R)  -  20  log  \F\.  (6) 

The  quantity  PF  is  often  called  the  propagation  factor,  while  PL  is  denoted  the  path  loss. 
The  VTRPE  code  outputs  either  PF  or  PL  values  on  a  user  specified  range/altitude  grid. 

The  remainder  of  this  report  is  divided  into  sections  which  cover  the  VTRPE  model 
physics  and  its  operation.  Section  2.  on  the  VTRPE  model  physics,  should  be  read  before 
attempting  to  use  the  model.  The  VTRPE  code  inputs  are  described  and  sample  input  run 
streams  are  provided. 


i  Accession  For 
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§2.0  Model  Physics 


This  section  briefly  describes  the  physics  and  numerical  methods  used  in  the  develop¬ 
ment  of  the  VTRPE  code,  including  the  split-step  Fourier  PE  (SSFPE)  algorithm.  Addi¬ 
tional  details  and  references  may  be  found  in  a  companion  report  by  Ryan.3  This  section 
is  designed  to  help  the  reader  better  understand  the  usage  of  various  VTRPE  model  inputs 
that  are  covered  in  later  sections. 

First,  the  electromagnetic  vector  wave  equations  are  derived  from  Maxwell’s  equations 
and  then  reduced  to  scalar  Helmholtz  form  via  an  earth-flattening  coordinate  transfor¬ 
mation.  Second,  the  PE  approximation  to  the  reduced  wave  equation  is  derived  and  its 
solution  via  the  SSFPE  algorithm  is  discussed.  Third,  the  calculation  of  the  pattern  propa¬ 
gation  factor  using  free-space  dipole  fields  is  reviewed.  Finally,  various  details  of  the  SSFPE 
algorithm  implementation,  including  starting  fields  and  range  step  sizes,  are  covered. 

2.1  Vector  Wave  Equation 


To  calculate  the  complex  pattern  propagation  factor  F.  Eq.  (4).  requires  knowledge  of 
the  macroscopic  electric  E  and  magnetic  H  radiation  fields.  These  fields  are  solutions  to 
vector  wave  equations  that  are  derived  from  Maxwell's  equations. 

In  a  source-free  region,  the  monochromatic  Maxwell’s  equations  (in  rationalized  inks 
units)  which  govern  the  propagation  of  electromagnetic  waves  are3-0 

V  x  E(r.cc’) 

V  x  H(r.u,') 

//0V  •  H(r.u.-) 

e0Y  •  f(r)E(r.u,’) 


=  H(r,w’). 

=  —  /u,’e05(r)E(r.  u.'). 
=  0. 

=  0. 


where  u;  is  the  radian  frequency  (u.’  =  2t r/)  and  the  implicit  time-dependence  exp (— /Mf) 
has  been  suppressed.  The  medium  electrical  properties  are  specified  via  the  free-space 
magnetic  permeability  //0  and  the  spatially  varying,  complex  relative  dielectric  constant 
f (r ).  For  microwave  frequencies  in  nonionized  media,  the  complex  dielectric  constant  may 
be  represented  in  the  form 


Mr) 


M(r)  +  > 


ofr) 
u-?eo  ' 


where  e0  is  the  vacuum  dielectric  constant,  a  is  the  medium  conductivity,  and  £]  is  the 
usual  permittivity  of  the  medium.  A  related  parameter  is  the  vacuum  wave  number  = 
u,’ %/// o e o  =  ^ /c-  "'here  c  is  the  speed  of  light. 

Using  standard  vector  identities.'  the  first  order  Maxwell’s  equations  (7)  may  be  re¬ 
placed  by  equivalent  second  order  vector  wave  equations  for  either  the  electric  or  magnetic 
field.  If  the  electric  field  E  is  eliminated  from  Eq.  (7).  then  the  magnetic  intensity  vector 
H(r)  satisfies: 

Y2H(  r)  4-  x  V  x  H(r)  +  £nc(r)H( r)  =  0  .  (S) 

Mr) 
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The  electric  field  is  then  determined  from  the  Maxwell  curl  relation 

E(r,u.-)  =  - - - V  x  H(r.u;). 

^e0c'(r) 

In  a  similar  fashion,  if  the  magnetic  field  H  is  eliminated  from  Eq.  (7).  the  electric  field 
vector  E(r)  may  be  shown  to  satisfy: 

V2E(r)  +  V  (E(r)  •  +  A-2c(r)E(r)  =  0.  (9) 

2.2  Scalar  Wave  Equation 


The  vector  wave  equations.  Eq.  (8)  or  (9),  may  be  converted  into  simpler  scalar  wave 
equations  if  the  transmitter  emits  radiation  that  is  linearly  polarized —  i.e..  the  electric 
field  vector  has  nonzero  components  lying  either  wholly  within  (vertical  polarization)  or 
perpendicular  to  (horizontal  polarization)  the  meridian  plane  containing  the  source  and 
observation  point.  This  is  the  case  for  many  types  of  radar  systems. 

Calculations  are  done  in  a  spherical,  earth-centered  coordinate  system  (r.B.o).  with 
respective  unit  vectors  (£r.£g.t0).  In  this  coordinate  system,  one  of  the  linearly  polarized 
field  components  is  in  the  azimuthal  direction  £p.  For  vertical  polarization,  the  nonzero 
magnetic  field  component  is  H(r)  =  H0( r)£0  which  satisfies 


£_d_l  d(rH0)  e  d  sin B  dHQ 
r  dr  s  dr  ^  r2  sin  6  SB  e  dB  + 


cot  0  dz 

~^7d0 


H0=  0. 


For  horizontal  polarization,  the  nonzero  electric  field  component  is  E(r)  =  E0(r)£0  and 
satisfies 

\d2{rE0{r))  1  d  .  dE0(r) 

r  dr2  +r2sinfl<^Sm  d9  + 


A’oc(r)  - 


•2  sin2  0 


Ea(  r)  =  0. 


The  above  equations  for  H0  or  Ea  are  then  converted  to  the  scalar  Helmholtz  form 

?r(.r.  c)  =  0.  (10) 


d2  d2 

+  77-9  +  A  (r.z) 


d.  r2  d~ 2 

by  means  of  an  exact  earth  flattening  transformation:8 


.r  =  ad. 

z  =  oln(l  +  h /a). 

where  a  is  the  effective  earth  radius,  and  h  =  r  —  a  is  the  local  altitude.  A  new  dependent 
variable  w  is  then  defined  by 

horizontal  polarization:  «'/,  =  Vr  sin  B  E0(  r). 

_  f+:/(2o)  ^Jg  sin(j -/a)  E0(x). 

/rf0(x).  if  x/a  <C  1  and  z/a  <C  1: 


(11) 


and  bv 


•  ,  ,  .  Vr  sin  6  TT 

vertical  polarization:  u\.  =  - 


n 


_  +3: /(2g)  y/asm(x/a) 


m(x.z) 


HM. 


n. 

- H0[x).  if  xjn  <C  1  and  z/a  <C  1. 


A  new  “effective"  wave  number  A  is  also  defined  bv 


horizontal  polarization: 
vertical  polarization: 


A’?  = 

I<1  =  him2 


3s  ec2{x /a) 
4  a2 

3  sec2(x/a) 
4a'2 


/  d2m  1 


cot(  r/a)  dm 

a7 

dm 


am 

d2m~i 

dz2 


adz 


where  m  is  the  modified  index  of  refraction 


m(.r,  c)  =  n(r) 


(12) 


(13) 


(14) 


2.3  Boundary  Conditions 


Maxwell's  equations  are  augmented  with  boundary  conditions  on  the  field  components. 
These  are  (1)  a  Sommerfeld  radiation-tvpe  boundary  condition  at  infinity9 


lim  r 

- - *OC 


(15) 


where  A  denotes  H0  or  E0.  and  (2)  continuity  of  the  tangential  electric  and  magnetic 
fields  at  the  earth's  surface  r  =  a  .  Continuity  of  tangential  field  components  is  achieved 
by  modeling  the  earth  as  a  locally  homogeneous  dielectric  with  finite  conductivity  and 
specifying  a  surface  boundary  condition.  This  surface  boundary  condition  is  implemented 
via  the  Leontovich  surface  impedance  boundary  condition10 


d{  rA) 


rdr 


=  -ZA 


where  the  quantity  Z  is  the  local,  flat-surface  impedance  defined  by 


-  1 


for  horizontal  polarization, 
for  vertical  polarization. 
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with  cs  the  complex  dielectric  constant  of  the  earth's  surface. 

The  Leontovich  impedance  boundary  condition  in  earth-flattened  coordinates  becomes 


3irhiT.z] 

a z 

for  horizontal  polarization,  and 
du'v(x.  z 


c=0 


—  +  >koy/6s  ~~  1  )  r/i(J’-0) 


^  I  c=0 
for  vertical  polarization. 


/  Jl  1  dm(x.z) 

\  2o  m(j.O)  dz 


+  ik o - 7 -  j  u\  (x.  0) 


:=0 


2.4  Pattern  Propagation  Factor 


If  the  transmitter  emits  only  linearly  polarized  radiation,  the  pattern  propagation  factor 
jF,  Eq.  (4).  only  involves  the  field  component  along  <0: 


i^v.r 


horizontal  polarization 


vertical  polarization 


where  E ^  and  are  the  o-component  of  the  free-space  electric  and  magnetic  Helds 
respectively.  To  calculate  the  pattern  propagation  factor  F.  the  free-space  fields  must  first 
be  specified.  By  convention,  the  pattern  propagation  factor  is  defined  with  respect  to  the 
free-space  field  of  a  unit-strength  point  dipole. 

For  horizontal  polarization,  the  reference  free-space  electric  field.  Emj.  is  chosen  to  be 
that  arising  from  a  vertically  oriented  magnetic  dipole  ( VMD).  For  vertical  polarization,  the 
reference  magnetic  field.  H€C/.  is  that  from  a  vertically  oriented  electric  dipole  (YED).  For 
both  the  YED  and  YMD  cases,  the  unit-strength  point  dipole  is  located  at  Tq  —  (r0.0.0) 
in  a  spherical  earth-centered  coordinate  system  (r.B.o)  with  the  dipole  moment  oriented 
along  the  polar  axis. 

Employing  the  dyadic  Green's  function,  the  o-components  of  the  dipole  fields  are  then 
given  by 1 1 


(1G) 


(IT) 


where  R  =  |r  —  r0|  =  y  r2  +  —  2rr0cos#  is  the  radius  vector  between  the  dipole  and 

the  field  observation  point,  and  ft  is  the  polar  angle. 


Employing  Eq.  (17)  and  (11),  the  propagation  factor  for  horizontal  polarization  F/(  is 
computed  as 


Fh(r)  = 


F0(r) 

Emd(r) 

4tt  \irh{r)\R2 
^O^V'o  ( r  sin  0)  ? 


1  +  (  krxR  ) 


Similarly,  using  Eq.  (1C)  and  (12).  the  propagation  factor  for  vertical  polarization  Fv  is 


Fr(  r)  = 


H0(r) 

Hed{r)  ' 

47t  |n(r)ir,.(r)|J?2 
^0U'’  (rsin#)2 


1  +  ( A'0f?) 


In  earth-flattened  coordinates,  the  propagation  factor  takes  the  following  limiting  form 
for  altitudes  and  ranges  typical  of  tropospheric  propagation: 


FhW 


F(.(x) 


4 *r2  x-$  in(x)l 

^O^/'o  \/l  +  (k0R)~2 

4tTv/j  \wh(  X;j  ,f  ( 

^0U'7,0  \/l  +  U'o>r)-2 
47r/?2  _ 3  |m(x)u\.(x)| 

_  .  o  _ _ _* _ __ 

V’  *  V1  +  (W2' 

4?r  v/x|m(x)trr(x)| 


~o )/J  «  1 


^0*’  \/l  +  ( A'o  >r )  2 


>  -  -o  )/.r  <  1 


2.5  Split-Step  PE 

The  propagation  factor  F  requires  the  field  O-component .  which  is  directly  related  to 
the  solution  i r  of  the  scalar  wave  equation  (10).  To  calculate  w.  the  elliptic  Helmholtz 
equation  is  first  approximated  by  a  parabolic  wave  equation  that  in  turn  is  solved  via 
formal  operator  methods.  Making  the  envelope  transformation 


u'(  x )  =  f?AVV(x). 

with  k'Q  the  reference  wave  number,  the  elliptic  wave  equation  is  replaced  by  the  generalized 
one-way  parabolic  wave  equation  (GPE) 


dv(x) 

d.r 


iQ{x)i'(x) . 


H 


The  pseudodifferential  GPE  operator  Q  is  defined  by 

<?(x)  =  \jd2/t ):2  +  A'2(x)  -  k0. 

Next,  the  "wide-angle"  split-operator  PE  approximation  to  Q  is  made  in  terms  of  ordi¬ 
nary  operators:12 


Q(x.z)*  A(z)  +  B(. r,c). 


A(:)  =  +  a2/a--2  -  «■ 

B(x.z)  =  K(x.= )  -  A-0  . 


0  • 


Then  a  formal  solution  of  Eq.  (IS)  can  be  expressed  in  exponential  operator  form  via  the 
Magnus  expansion13 

+?  A. AA  +  B) 


\'(X  +  A,..; 


where  B  =  B(.r  +  jA^.r). 

The  Trotter  product  formula1"1  is  then  used  to  symmetrically  factor  the  Magnus  expan¬ 
sion  into  the  product  of  simpler  operators.  L  ;  and  UB.  where 

r.,  =  r.,(  = 

l-B=UBiA,)  =  ,+iA*B. 
to  yield  the  split-step  PE  (SSPE)  algorithm: 


v(x  +  Ax.z)  =  UA(Ar/2)UB(AI)UA(Ar/2)v(x.z). 

=  f +?Ax.4/2  (+?Aj.B  (+iArA/2  L^x  ~^ 


(19) 


The  PE  solution  at  range  x  +  Ax.  v{x  +  Ax.  c ).  is  obtained  from  the  known  field  at  range 
x .  v{  x. :).  by  means  of  the  SSPE  algorithm.  Physically,  the  SSPE  algorithm  amounts  to 

(1)  a  half-step  of  free-space  propagation.  l\(Ax/ 2). 

(2)  a  phase  correction.  UB( Ax).  to  account  for  refractive  effects,  and 

(3)  another  half-step  of  free-space  propagation.  U  a{At/2). 


2.6  Split-Step  Fourier  PE  Algorithm 


The  presence  of  the  differential  operator  .4( c )  in  the  SSPE  exponent.  Eq.  (19).  is  dealt 
with  by  transforming  to  a  basis  in  which  .4(c)  is  diagonal.  One  such  basis  is  the  Fourier 
basis.  The  c-space  Fourier  transform  T  is  defined  as 


$  ( x .  p ) 


^■[v(jr.c)]  =  f  vix.z) f  !p:  d:. 

J—O C 
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with  the  corresponding  inverse  transform  T  1  being  defined  by 


t'(i\z)  =  T  l[H '{r.p)}  =  -^-f  ^(.r.p)(+,pz  dp. 

1/  —  'DC 

The  transform  variable  p  may  be  associated  with  a  vertical  wave  number  via  p  —  k0  sin  8. 
where  8  is  the  propagation  angle  with  respect  to  the  horizontal. 

The  exponential  operator  f is  evaluated  using  Fourier  transform  methods  as 


c.iA,*4(c)v(;r  =)  =jr-l 


f  “*Ar  (*0  -  \Jkl~P 2) 


(20) 


Employing  Eq.  (20)  in  Eq.  (19)  yields  the  split-step  Fourier  PE  (SSFPE)  algorithm: 

I -  a  B  ,  \  ~'Ax  (ko  -  \] H  ~  P2]  r 

V(*  +  Ar;)  =  (+,^JJF“1  {(  V  /  JP[t>(.r.c)] 


(21) 


Providing  the  initial  starting  field  e(0,r)  is  known,  the  SSFPE  algorithm  can  efficiently 
advance  the  solution  in  range. 

2.7  PE  Starting  Fields 

The  split-step  PE  method  must  be  initialized  with  a  starting  field  distribution  xl’{ 0.  ;). 
This  is  accomplished  in  the  VTRPE  model  by  using  the  duality  of  the  e-space  antenna 
aperture  field  distribution  -4(e)  and  the  corresponding  p-space  radiation  pattern  func¬ 
tion  f(p).l° 

In  free  space,  the  antenna  radiation  pattern  fimetion  f(p)  and  the  antenna  aperture 
field  distribution  -4(e)  are  a  Fourier  transform  pair: 


f(p)=  r  dz . 

J  —  DC 


where  the  p-space  parameters  are  defined  by 

p  =  k-Q  sin  8. 

p„  =  k()  sin#,)-  (22) 

In  Eq.  (22).  8  is  the  elevation  angle  measured  with  respect  to  the  horizontal  (8  >  0  is 
up)  and  8 0  is  the  antenna  mainlobe  vertical  pointing  direction.  The  antenna  pattern 
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function  f{p)  is  used  to  construct  even.  \k£(0 .p).  and  odd.  ^o(0,p).  symmetry  p- space 
starting  fields  in  the  form 


*e(0  .p)  =  f(p)  +  f(-p)e+iP:o. 


(23) 


and 


*o(0.p)  =  f(p)(-ip:o  -  f(~p)e+ip:°. 


(24) 


where  the  Fourier  shift  theorem  has  been  used  to  properly  account  for  a  nonzero  source  al¬ 
titude,  c0.  Given  the  above  even/odd  symmetry  p-space  field,  the  corresponding  even/odd 
symmetry  c -space  field  is  obtained  by  taking  the  inverse  cosine  or  sine  transform  of  Eq.  (23 ) 
or  Eq.  (24).  respectively: 


ve(0.=)  =  T-l[<be(0.p)}, 
Vo(0.r)  =  jF-1[^o(0  ,p)}. 


In  practice,  the  infinite  Fourier  transforms  in  Eq.  (25)  are  replaced  by  finite  discrete 
sine  or  cosine  transforms  over  the  finite  interval  (0.  Zniax  ).16  These  discrete  transforms  are 
in  turn  evaluated  numerically  using  real- valued  fast  Fourier  transform  (FFT)  methods.1' 
Use  of  FFTs  requires  proper  selection  of  the  transform  size  -V  and  the  sampling  interval 
Zniax  to  ensure  satisfaction  of  the  Nyquist  theorem  and  thereby  avoid  aliasing  problems.18 
If  p  =  k'Q  sin0max  is  the  maximum  vertical  wave  number  in  the  FFT  ;> space  transform 
grid,  and  X  is  the  discrete  cosine/sine  transform  size,  then  the  relation  between  Zmax  and 
p  is  given  bv 

'max  ° 

Zm  a  v  P  —  7T  J\  . 
max'max 


Equivalently,  the  specification  of  the  FFT  :r-grid  size.  6,.  and  the  corresponding  p-grid 
size.  6p,  is  provided  by 


6~  =  7 r/p 

'  ,nax 

^ p  =  ^/^max  • 


(26) 


Two  constraints  need  to  be  placed  on  the  FFT  p- space  grid  to  ensure  accurate  results 
using  the  SSFPE  method.  These  constraints  depend  on  the  initial  p- space  field  spectrum, 
and  on  the  desired  vertical  resolution  in  c-space.  First,  the  value  of  p  must  be  suffi- 
ciently  large  to  contain  all  significant  energy  in  the  starting  field  'P(O.p).  If  's  the 

elevation  angle  corresponding  to  a  specified  sidelobe  level  in  the  antenna  radiation  pattern 
(e.g..  50  dB  down  relative  to  the  mainlobe  peak)  and  pSJ  =  k0  sin &SLL-  an  estimate 
for  p  is 

max 

Pmax-*W  +  W-  (27) 

Typically,  only  sidelobes  above  a  certain  threshold  are  included  in  the  VTRPE  code  to 
minimize  the  size  of  p  .  Before  transforming  back  from  p-space  to  c-space.  the  initial 
field  ^(O.p)  is  low-pass  filtered  to  remove  energy  above  the  Nyquist  point  Pmax  thereby 
avoiding  FFT  aliasing  problems.  Since  p  is  related  to  the  vertical  spacing  d.  via 
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Eq.  (26).  a  further  constraint  on  its  size  is  provided  by  the  smallest  vertical  scale  length  lr 
needed  to  resolve  the  refractivity  field  N{z): 


>*/lv 


(28) 


The  maximum  vertical  wavenumber  p  is  a  function  of  the  antenna  pattern  sidelobes. 

max 

and  the  vertical  tilt  angle  #0. 

A  second  contraint  is  placed  upon  the  p-space  mesh  size  dp.  To  adequately  resolve  the 
surface  image  interference  lobes  in  the  starting  field  'I'(O.p),  the  p- space  resolution  tp  must 
satisfv 


bp  z 


r—< 

max 


This,  in  turn,  constrains  the  minimum  allowable  size  of  Zma,.  to  be 


7  >  o~ 

max  ^  "~0 


2.8  Analytic  Antenna  Patterns 


The  SSFPE  method  is  capable  of  modeling  the  radiation  emitted  by  directional  an¬ 
tennas — provided  that  the  complex  antenna  radiation  pattern  f(p)  is  known.  Often  this 
information  is  not  available  for  specific  radar  systems,  so  simplified  generic  antenna  pat¬ 
terns  are  used  instead.  These  generic  patterns  display  some  of  the  features  of  real  antennas 
while  retaining  fairly  simple  analytical  forms.  Each  of  the  analytic  antenna  patterns  is 
specified  in  terms  of  a  normalized  p-space  steering  parameter  t  defined  by 

sin  9  —  sin  90  P  ~  P0 
sin  (?9bw)  p i. 

where 

PL  =  *osin(-06u.). 

and  9f)W  is  the  half-power  beamwidth  (HPBW).  In  addition  to  the  HPBW.  each  antenna 
pattern  has  sidelobes  specified  at  a  certain  power  level  relative  to  the  mainlobe  peak. 

The  following  analytic  antenna  radiation  patterns  are  useful  approximations  to  real 
antenna  radiation  patterns  and  are  implemented  in  the  VTRPE  model: 

2.8.1  sin(x)/x 


The  uniformly  illuminated  aperture  corresponds  to  a  radiation  pattern  having  the  func¬ 
tional  form 

„  sin(of) 

f(t)= — - — -.  n  =  1.391C. 

J  at 

The  scale  parameter  a  is  determined  by  solving  the  nonlinear  equation 


sin  o 


v/2  * 


This  pattern  has  the  narrowest  mainlobe  width,  at  the  expense  of  high  sidelobe  levels. 
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2.8.2  Gaussian 


The  Gaussian  antenna  radiation  pattern  has  the  functional  form 

o  hi  2 

f(t)  =  (~at.  where  a2  =  — . 

This  pattern  is  the  optimal  compromise  between  sidelobe  level  and  mainlobe  beamwidth. 

2.8.3  Compound 

The  compound  radiation  pattern  is  a  one-parameter  pattern  that  is  a  blend  of  the 
uniformly  illuminated  aperture  and  the  cosine- squared  aperture  distribution,  having  the 
functional  form 


m  = 


2  sin(af) 
1  4-  C b  at 


1  -Q  1 

2  1  -  (at/-)2 


0  <  Cb  <  1. 


The  parameter  a  controls  the  mainlobe  width  and  is  determined  by  solving  the  nonlinear 
equation 

i  +  O 

v/2  ' 

The  uniform  aperture  corresponds  to  C(,  =  1.  while  the  cosine-squared  aperture  corre¬ 
sponds  to  Of,  =  0.  The  cosine-squared  aperture  is  often  used  as  an  approximation  to  the 
radiation  pattern  from  a  parabolic  reflector  antenna. 


sm  a 


20,  + 


i-o 


1  -  (a In)'2  j 


2.8.4  Hansen 


Another  single  parameter  pattern  is  obtained  from  the  circular  aperture  distributions 
analyzed  by  Hansen.19  This  pattern  has  the  functional  form 


where 


m  = 


— 
',(«)  i 


V 


?i(\/0 

\/M 


for  .r  >  0 
for  x  <  0. 


x  =  H 


2 


(at)2. 


and  j  j  and  tj  are  the  first-order  spherical  and  modified  spherical  Bessel  functions  respec¬ 
tively.  The  quantity  a.  which  determines  the  3-dB  point  in  the  pattern,  is  found  by  solving 
the  transcendental  equation 


The  first  sidelobe  level,  SLL.  in  the  radiation  pattern  can  be  shown  to  have  a  value  of 


SLL  =  30.84  +  20  log 


3m  (H) 

H 


dB. 


down  from  the  peak,  and  is  located  at  p  =  p 


SLL 


Ps Li.  ~Pq  =  ^-yJl  +  (o.7G3/H)2. 

The  parameter  H  allows  a  trade-off  between  low  sidelobe  level  and  mainlobe  beamwidth. 


2.9  SSFPE  Range  Step  Size 


Each  application  of  the  SSFPE  algorithm.  Eq.  (19).  leads  to  a  local  truncation  error 
in  the  solution  r  that  is  proportional  to  the  cube  of  the  PE  range  step  A_,..  Since  many 
range  steps  are  typically  taken,  these  local  errors  can  accumulate  and  produce  unacceptable 
errors  in  the  solution.  To  prevent  this  from  happening,  the  VTRPE  code  performs  a  global 
error  estimate  based  upon  a  detailed  analysis  of  the  local  SSFPE  errors.  To  bound  the 
total  global  error  in  the  solution  the  PE  range  step  A..  is  chosen  so  that 


2.„,.  II*"’-1 


<  24 

where  erei  is  a  specified  error  tolerance,  and 

r^max 


><"c| 


2*0  Jo 


V 


—  -  k2V 

dz2  0 


rdz. 


II*(3V|| 


48A-, 


(  max 

wa,"Y  M’ + a2r 

dr 

2  j  d2V  d2\r\2' 

[Jo 

°\dz)  K  1  +  dz 2 

dz 

4  dz2  dz2 

dz 


+ 


l  d2V 


dz* 


;=0 


dz 


Here  F(~)  =  m2{z)  —  1  is  the  "effective  potential"  related  to  the  modified  index  of  refraction 


m. 


The  SSFPE  range  step  A^  is  then  set  by 


24e„ 


A,  = 


f  ^max 

/  ^ 

J  o 


91  _  ,2,. 

3-2  ‘o' 


V  I 


f2  idv\2.  |2  d2v  dr 

kn  (  ~x~  I  (t'|  +  q  2 


t  *••(£) 


!  d2v  d2\v\2 
4  dz 2  dz 2 


*  +  T 


d2V  d\i 


]  V  t  u\ l ' 

~d^~dz 


:=0 


As  the  VTRPE  code  advances  the  field,  the  local  error  budget  is  monitored  and  the  range 
step-size  A_,.  is  dynamically  adjusted  to  keep  the  local  error  below  a  preset  threshold. 
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2.10  Reflectionless  Absorber 


The  appropriate  boundary  condition  to  be  satisfied  as  c  — +  oc  by  the  PE  field  r(.r.  c)  is 
the  Sonimerfeld  outgoing  wave  radiation  condition.  Ecj.  (15).  Since  the  split-step  Fourier 
algorithm  employs  finite  Fourier  transforms,  the  implementation  of  a  radiation-type  bound¬ 
ary  condition  is  quite  complicated.  This  follows  from  the  fact  that  truncation  of  the  infinite 
c- domain  down  to  a  finite  interval  in  the  Fourier  transform  leads  to  the  introduction  of 
spurious  discrete  standing  wave  solutions  in  the  vertical.  In  effect,  the  terminal  impedance 
at  the  end  of  the  transform  grid  is  not  properly  “matched"  to  the  radiation  boundary 
condition. 

To  circumvent  this  problem  and  attenuate  the  spurious  standing  wave  solutions  intro¬ 
duced  by  the  finite  Fourier  transforms,  a  complex  absorber  potential  V(I^s(c)  is  added  to 
the  split-step  B  operator: 


B(j-.z)  =>  B(.r.:)  +  ~Vab ,s.(c). 

The  specific  form  chosen  is 

Vabs(  =  )  =  1  o  seclr [( c  -  Zinax  )/«•„]. 

where  the  parameters  {Vo-11'*)}  are  determined  parametrically  by  minimizing  transmission 
and  reflection  coefficients  from  the  sponge  region.20  The  actual  "physical"  propagation 
region  in  the  VTRPE  model  thus  extends  to  an  altitude  -  =  /fmax  <  Z)llax  to  allow  room 
for  the  absorber. 


§3.0  Usage 


This  section  describes  how  to  access  and  run  the  VTRPE  code.  It  assumes  that  the  VTRPE 
code  has  already  been  installed  on  either  a  DOS-based  PC  or  UNIX-based  system,  and  that 
the  VTRPE  executable  is  in  the  current  working  directory.  Optionally,  an  environmental 
path  variable  may  be  set  to  the  directory  containing  the  VTRPE  executable  (this  is  the 
preferred  method,  since  it  allows  sub-directory  organization  of  data  files). 


3.1  Notation 

The  following  notation  conventions  are  used  in  this  report. 

3.1.1  Typographic  Conventions 

The  following  typographic  conventions  are  used  in  later  sections  of  this  report  to  describe 
various  input  formats. 

{ }  Braces  are  used  to  group  items  together. 

()  An  item  within  angle  brackets  (including  the  brackets)  is  a  parameter  which 
needs  to  be  replaced  by  user- supplied  text. 

[  ]  Anything  within  square  brackets  is  optional:  the  brackets  are  not  part  of  the 
bem. 

The  ellipsis  means  that  the  previously  defined  pattern  can  be  repeated  several 
times.  The  ellipsis  represents  the  position  at  which  text  or  items  have  been 
omitted. 

text  Input  text  to  be  entered,  for  example  program  names  or  input  keywords,  is 
shown  in  a  monospaced  text  font. 

3.1.2  Namelist  Input 

The  VTRPE  inputs  are  invoked  via  ASCII  files  in  FORTRAN  namelist  format.  This 
type  of  file  consists  of  plain  text  records  up  to  80  characters  in  length  (i.e.."card  images"  ) 
and  may  be  prepared  using  a  text  editor  or  word  processor.  If  a  word  processor  such  as 
WordStar  or  WordPerfect  is  used,  be  sure  to  disable  any  “word  wrap  or  special  output 
formatting  options. 

FORTRAN  namelist  input  is  a  powerful,  free-format  type  of  input,  consisting  of  ordered 
strings  of  “items"  or  symbolic  variable  names,  and  their  associated  "values."  An  item  is  a 
mnemonic  name  (e.g..  MEGAHZ  for  the  transmitter  frequency  in  MHz)  that  is  associated  with 
a  particular  FORTRAN  variable  scalar,  array  name,  or  array  element.  Associated  with 
each  item  in  the  namelist  input  stream  is  a  value,  that  is  a  FORTR AN-type  INTEGER. 
REAL.  DOUBLE  PRECISION.  COMPLEX.  LOGICAL,  or  CHARACTER  constant  that 
is  appropriate  for  the  associated  item. 

Different  types  of  namelists  are  used  to  specify  different  categories  of  inputs  to  VTRPE. 
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Each  namelist  input  field  has  the  following  general  form: 

kzident 

namelist  data  fields 

&END 

The  start  of  a  namelist  input  stream  is  signaled  by  a  unique  identifier  string  of  the  form 
"kident'  .  where  ident  is  the  name  of  a  particular  namelist  (e.g..  &RPE).  while  the  end  of 
the  input  stream  is  signaled  by  the  unique  terminator  string  ''&END".  In  between  the  initial 
and  terminating  identifier  strings  is  the  body  of  the  namelist,  consisting  of  free-format  data 
fields.  The  following  rules  apply  to  the  namelist  identifier  strings: 

•  The  initial  identifier  string  kidcvt  must  start  in  column  two  of  the  input  record,  be 
capitalized,  contain  no  embedded  blanks,  and  be  contained  within  a  single  record. 

•  At  least  one  blank  space  must  separate  the  starting  identifier  string  from  the  body 
of  the  namelist. 

•  The  terminating  identifier  string  is  always  &END. 

Otherwise,  white  space  (i.e..  blanks)  may  appear  anywhere  in  the  list. 

Each  data  field  or  input  record  in  the  namelist  has  the  general  format: 

item  j  =  value  j  [<  sep  >  [?/em2]  =  value^}--- 

where  item  is  the  symbolic  name  of  a  scalar,  array,  array  element,  or  character  string:  value 
is  a  FORTRAN-type  INTEGER.  REAL.  DOUBLE  PRECISION.  COMPLEX.  LOGICAL, 
or  CHARACTER  constant  that  is  appropriate  for  the  associated  item:  and  <  sep  >  is  a 
comma,  tab,  or  space  separator.  The  actual  symbolic  names  and  associated  data  types  for 
various  namelist  items  are  covered  later. 

The  following  rules  apply  to  namelist  data  field  records: 

•  An  item  specified  on  the  left  side  of  the  equal  sign  cannot  contain  spaces  or  tabs 
except  within  the  parentheses  of  a  subscript.  Each  item  must  be  contained  in  a 
single  record.  The  spelling  and  case  of  the  item  name  is  critical! 

•  CHARACTER  string  constants  are  enclosed  in  single  quotes. 

•  COMPLEX  constants  are  specified  in  the  format  item  —  (real.  imag)  by  enclosing 
the  real  and  imaginary  parts  of  the  complex  quantity  in  parenthesis,  separated  by 
a  comma. 

•  A  null  value  is  specified  by  two  consecutive  commas,  by  an  initial  comma  or  by 
a  trailing  comma.  A  null  value  indicates  that  the  corresponding  namelist  array 
element  is  unchanged.  A  null  value  can  represent  an  entire  complex  element  but 
cannot  be  used  for  only  one  part  of  a  complex  constant . 

•  Multiple  items  may  appear  on  a  single  line. 

•  The  equal  sign  in  a  value  assignment  can  be  delimited  by  spaces. 

3.2  Program  Invocation 

To  invoke  the  VTRPE  model,  simply  enter  VTRPE  at  the  command  line  prompt: 
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>  VTRPE 

where  ">"  denotes  the  operating  system  command  line  prompt.  After  the  model  is  loaded 
into  memory,  a  logo  will  appear  and  the  user  will  be  prompted  for  the  names  of  the  input 
and  output  files  to  be  used  in  the  mil.  The  content  and  format  of  the  various  input  and 
output  files  are  described  in  subsequent  sections.  Optionally,  the  run  input  and/or  output 
filenames  may  be  specified  at  program  invocation  via  command  line  arguments  that  follow 
the  program  name: 

>  VTRPE  input  _f  ih  apt  c[.  print  .f  il(  spt  c ]  [.plot-fill  spt  c] 

These  run-time  arguments  consist  of  valid  filenames  (with  an  optional  path)  in  the  order 
indicated.  Items  inside  the  square  brackets  are  optional.  Each  fill. spec  denotes  the  name 
of  a  file  to  be  processed,  which  can  include  a  full  or  partial  path  specification.  Not  all 
filenames  need  be  specified  on  the  command  line;  those  that  are  not  will  be  prompted  by 
the  program. 

VTRPE  inputs  are  invoked  via  a  formatted  IXPUT-file  denoted  by  input./ ill  apif  .  The 
contents  and  format  of  the  IXPUT-file  are  described  later  in  this  document.  In  addition 
to  the  IXPUT-file.  there  are  two  optional  output  files:  (1)  the  PRIXT-file.  denoted  by 
print. fill  spic.  which  is  an  ASCII-formatted  file  containing  diagnostic  and/or  debug  out¬ 
put.  and  (2)  the  PLOT-file.  denoted  by  plot. f  ih spec,  which  may  be  either  ASCII  or  binary 
and  contains  a  compacted  range  altitude  matrix  of  loss  values  that  may  be  postprocessed 
and  input  into  a  user  specified  plotting  program.  The  file  formats  of  the  PRIXT  and  PLOT 
files  are  covered  later. 

In  addition  to  the  above-mentioned  files,  the  VTRPE  code  will  automatically  allocate 
temporary  “scratch"  files,  which  are  deallocated  upon  successful  run  termination.  If  any 
errors  have  occurred  during  program  execution,  then  a  file  having  the  same  filename  as 
input. fill  spic  but  with  the  extension  “.err"  will  be  created  in  the  current  working  direc¬ 
tory.  Consult  this  *  .err  file  for  more  details  of  any  errors  that  may  have  occurred. 

As  the  VTRPE  program  is  run.  many  different  input/output  files  are  created,  and  it  often 
proves  beneficial  to  create  separate  subdirectories  for  different  projects.  Any  valid  filename 
may  be  used  for  the  input/output  filenames.  However,  to  facilitate  file  management  and 
prevent  confusion,  a  uniform  file-naming  convention  is  strongly  recommended  for  the  input 
and  output  filenames.  One  possible  choice  is  to  use  a  unique,  descriptive  run  filename  and 
then  vary  the  file  extension  depending  upon  the  type  of  file.  For  example,  the  following 
default  file  extensions  are  suggested: 

jobnawi .inp  =>  IXPUT-file 
jobnami  .prt  =>  PRIXT-file 
jobn ami  .pit  PLOT-file 

3.3  Input  Files 

Inputs  to  the  VTRPE  model  are  via  ASCII  formatted-type  files  in  FORTRAX  namelist 
format.  Three  types  of  namelist  input  streams  are  processed  by  VTRPE  : 

(1)  RPE. 


(2)  SECTOR. 

(3)  GROUND. 


The  general  layout  of  the  VTRPE  input  stream  is  shown  in  figure  1. 

&RPE 

transmitter  characteristics  and  run-time  options 

&END 

&SECT0R 

:  refractivity  profile  data  sector  1 

&END 

&SECT0R 

:  refractivity  profile  data  sector  2 

&END 

&SECT0R 

refractivity  profile  data  sector  N 

&END 

&GR0UND 

:  surface  elevation  and  dielectric  properties  (optional) 

&END 


Figure  1.  VTRPE  input  stream. 

The  first  namelist.  RPE,  contains  information  relating  to  the  transmitter  characteristics 
(e.g..  frequency,  beamwidth,  or  polarization)  and  general  run-time  parameters  (e.g..  output 
range  stepsize).  The  RPE  namelist  should  occur  first  in  the  input  file.  Section  4  describes 
the  variables  in  the  RPE  namelist. 

Following  the  RPE  namelist  are  one  or  more  SECTOR-type  namelists  which  contain  en¬ 
vironmental  refractivity  profile  data.  For  a  range-independent  run.  only  a  single  SECTOR- 
type  namelist  is  required.  For  a  range- dependent  run.  a  separate  SECTOR  namelist  is 
required  for  each  distinct  environmental  refractivity  profile.  The  refractivity  profiles  are 
grouped  into  data  sectors  that  are  to  be  arranged  in  monotonically  increasing  great  circle 
range  from  the  transmitter.  Section  5  describes  the  variables  in  the  SECTOR  namelist. 

Finally,  if  variable  terrain  or  finite  surface  conductivity  conditions  are  to  be  modeled, 
then  the  surface  elevation  profile  and  corresponding  surface  dielectric  properties  must  be 
input  via  the  GROUND-type  namelist.  In  this  case,  the  variable  terrain  control  flag  IVTPE 
must  be  set  in  the  RPE  namelist  to  flag  the  run  as  a  variable  terrain  run.  Various  formats 
are  available  to  input  the  surface  dielectric  properties  and  are  discussed  in  section  6. 
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§4.0  General  Inputs 


This  section  describes  the  variables  in  the  RPE  namelist.  These  variables  are  used  to 
set  initial  VTRPE  run  parameters  that  are  not  functions  of  the  environment.  These  include 

•  transmitter  height  and  antenna  characteristics. 

•  specification  of  the  model  output  altitude/range  grid. 

•  range/altitude  units  flags. 

•  SSFPE  range  step  size  and  FFT  size,  and 

•  type  of  VTRPE  model  output. 

4.1  Output  Range/ Altitude  Grid 

The  VTRPE  code  generates  output  (PF  or  PL)  on  a  user  specified  range/altitude  grid 
of  Xr  range  points  and  X.  altitude  points. 

The  vertical  grid:  {cr(/).t  =  1.2 . .V. }.  is  constant  with  range  and  defined  by  the 

entries  in  the  ZR  input  variable  described  below.  At  each  output  range,  the  VTRPE  code 
interpolates  the  complex  field  from  the  equi-spaced  FFT  c -space  grid:  -  id..  i  = 

0.1 . .V}  to  the  user  specified  receiver  grid 

The  Xr  =  NRANGE  output  range  points.  {r;  i  =  1.2 . A,-}-  correspond  to  the  nearest 

VTRPE  range  step  to  the  equi-spaced  range  grid  { }  defined  by 

+  (?  ~  1  i  =  1.  2 . Ar. 

where  /?0  =  RFIRST  is  the  initial  output  range  and  =  RXGIXC  is  the  output  range  in¬ 
crement.  Note  that  the  actual  output  ranges  {r?  }  may  not  correspond  to  the  user  specified 
ranges  {/?,}  due  to  the  fact  that  the  VTRPE  code  man  use  a  variable  range  step  size.  The 
user  must  specify  at  least  two  items  in  the  triplet  { .  Xr.  A#}:  the  remaining  component 
of  the  triplet  is  computed  internally  by  the  model. 

For  a  range-independent  run,  (i.e..  only  one  input  refract ivity  profile,  no  terrain  eleva¬ 
tion  variability,  and  simplified  surface  boundary  conditions)  the  following  RPE  namelist 
variables  are  the  minimum  set  needed  to  specify  the  run: 

•  transmitter  antenna  mainlobe  vertical  beamwidth.BMWIDTH 

•  transmitter  frequency.  MEGAHZ 

•  at  least  two  of  NRANGE,  RFIRST,  RLAST 

•  specification  of  the  output  range  increment.  RNGINC 

•  specification  of  the  output  altitudes.  ZR 

•  specification  of  the  transmitter  height.  ZTRANS 

All  other  RPE  namelist  variables  will  take  on  their  default  values.  To  change  a  particular 
variable,  merely  specify  it  in  the  namelist  input  stream. 

The  list  on  the  following  pages  describes  the  RPE  namelist  input  variables.  Mandatory 
input  parameters  are  indicated,  with  all  other  input  variables  optional.  If  a  nonmandatory 
input  variable  is  not  specified  in  the  RPE  namelist,  then  the  internal  default  value  will  be 
used.  Note  that  some  input  variables  that  normally  are  optional,  become  mandatory  if 
certain  options  are  specified. 
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4.2  RPE  Namelist  Variable  Descriptions 

ALPHA  Atmospheric  volume  attenuation  or(c) 

Total  frequency  dependent  atmospheric  volume  microwave  attenuation  vs.  altitude 
o,.(c).  Arising,  for  example,  from  absorption  by  molecular  gases,  water  vapor,  or 
liquid  water  (rain).  Stored  as  FORTRAN'  2-D  array  ALPHA(j.  /)  that  is  input  as 
ordered  pairs:  {^t.  a,},  i  <  50.  where  ALPHA(l.Z)  =  +  .  ALPHA(2./)  =  a-,  i  <  50 
where  +•  is  the  /th  altitude  and  o  ■  =  o(  zt )  is  the  corresponding  attenuation 
(dB/km)  at  the  altitude  The  input  attenuation  values  {c^.o,}  must  start  at 
the  surface  and  be  ordered  in  terms  of  increasing  altitude.  The  attenuation  values 
should  correspond  to  the  actual  run  frequency  as  specified  by  MEGAHZ.  The  attenua¬ 
tion  profile  a(c)  is  modeled  by  piecewise  constant  layers,  and  is  range- independent. 
The  altitudes  {;■}  need  not  correspond  to  the  input  refract ivity  profile  altitudes. 
If  volume  attenuation  is  desired,  then  also  specify  IVOLUME  variable  in  addition. 

UNITS:  attenuation  o;  in  dB/km:  altitude  +•  in  units  as  specified  by  IZUNIT 
variable 

DEFAULT:  ALPHA ()  =  0 

BMWIDTH  Transmitter  vertical  beam-width 

Free-space  1/2-power  (i.e..  angular-  width  to  3-dB  power  points)  transmitter  ver¬ 
tical  beamwidth.  Used  in  conjunction  with  mainlobe  tilt  parameter  BMELEV  to 
compute  the  initial  PE  starting  field.  The  particular  type  of  antenna  radiation 
pattern  to  be  used  in  generating  the  starting  field  is  specified  by  variable  ISOURCE. 
BMWIDTH  <  180.  use  free-space  pattern  having  specified  1/2-power  beamwidth 
BMWIDTH  =  ISO.  use  omnidirectional  point  source  pattern 
RANGE:  0  <  BMWIDTH  <  180 
DEFAULT:  Note:  mandatory  input 
UNITS:  deg 

BMELEV  Trans  mitter  vertical  mainlobe  beam  tilt 

Transmitter  free-space  antenna  pattern  mainlobe  vertical  tilt  angle  (with  respect 
to  the  horizontal). 

BMELEV  <  0.  transmitter  mainlobe  is  steered  down 
BMELEV  =  0.  transmitter  mainlobe  is  steered  horizontally 
BMELEV  >  0.  transmitter  mainlobe  is  steered  up 
RANGE:  -90  <  BMELEV  <  +90 
DEFAULT:  BMELEV  =  0 
UNITS:  deg 

CBEAM  Compound  aperture  parameter 

Controls  selection  of  the  PE  starter  aperture  field  -4(;).  which  varies  continuously 
from  a  uniform  aperture  (Q,  =  1)  to  a  cosine-squared  type  aperture  (Cf,  =  0) 


21 


depending  upon  the  value  of  Q,  =  CBEAM.  The  uniform  aperture  field  (CBEAM  =  1 ) 
corresponds  to  a  sin(  ,r)/.r  type  radiation  pattern,  while  the  cosine-squared  aperture 
field  (CBEAM  =  0)  is  a  common  approximation  for  parabolic  reflector  antennas.  If 
CBEAM  is  specified,  then  the  user  must  also  specify  I  SOURCE  =  1. 

CBEAM  =  0,  cosine  squared  aperture  field 
0  <  CBEAM  <  1.  compound  type  aperture  field 
CBEAM  =  1.  uniform  aperture  field 
RANGE:  0  <  CBEAM  <  1 
DEFAULT:  0 
UNITS:  none 

DBPASS  PE  wave  number  ■pass-band 

Controls  the  PE  wave  number  (p- space)  energy  spectrum  and  starter  field  pass- 
band  cutoffs.  The  initial  PE  fields  vertical  wave  number  power  spectrum  is  low- 
pass  filtered  for  levels  >  DBPASS  down  from  spectrum  peak.  DBPASS  is  also  used 
to  low-pass  filter  (i.e..  band  limit)  the  PE  field  at  each  range  step  and  controls 
variable  range  step  size.  May  be  used  in  lieu  of  #niax  =  THBMAX  to  limit  maximum 
PE  wave  number  pmax.  WARNING:  Do  not  specify  this  parameter  unless  sure  of 
results. 

DBPASS  =  0.  truncate  p-space  PE  field  based  upon  value  of  THBMAX  at  pmax  = 
*'0  sin  ^max 

DBPASS  >  0.  truncate  p-space  field  at  specified  sidelobe  level 
RANGE:  DBPASS  >  0 
DEFAULT:  DBPASS  =  50 
UNITS:  dB 

DRMAX  Maximum.  PE  range  step  size 

Used  in  variable  step  size  run  to  limit  the  maximum  SSFPE  range  step  size  that 
will  be  used.  Do  not  specify  larger  than  RNGINC  or  missed  output  will  result. 

RANGE:  0  <  DRMAX  <  RNGINC 
DEFAULT:  DRMAX  =  1000 
UNITS:  in 

DRMIN  Minimum  PE  range  step  size 

Used  in  variable  step  size  run  to  limit  the  minimum  PE  range  step  that  will  be 
used. 

RANGE:  DRMIN  >  0 
DEFAULT:  DRMIN  =  75 
UNITS:  meters 
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EPSURF  Complex  surface  dielectric,  constant 

Complex  effective  surface  dielectric  constant  at  the  transmitter  cs  (cs  =  cj  +  ?c2). 
If  EPSURF  =  0.  then  perfect  conductor  surface  boundary  conditions  are  used  in  the 
calculation:  if  nonzero,  then  finite  conductivity  surface  boundary  conditions  are 
used.  For  range-dependent  problems,  the  surface  dielectric  constant  cs.(/  )  is  input 
via  the  GROUND  input  stream.  The  complex  input  value  es  is  specified  as  an 
ordered  pair:  EPSURF  =  (c  j,52)-  where  fj  is  the  surface  permittivity,  and  f2  is  the 
relative  loss  function  (note  that  =  2  >0). 

EPSURF  =  (0..0.),  flags  perfect  conductor  surface 
EPSURF  /  (0..0.).  flags  finite  conductor  surface 
RANGE:  cj  >  1.  c2  >  0 
DEFAULT:  EPSURF  =  (0..0.) 

UNITS:  none 


FPLOT  PLOT-filename 

Character  string,  enclosed  by  apostrophes,  for  output  PLOT-filename.  If  non¬ 
blank.  then  the  specified  filename  is  used  for  PLOT-filename:  otherwise,  prompt 
for  filename  OR  use  command-line  argument  . 

DEFAULT:  FPLOT  =  blank 
UNITS:  none 


FPRINT  PRINT -filename 

Character  string,  enclosed  by  apostrophes,  for  output  PRINT-filename.  If  non¬ 
blank.  then  the  specified  filename  is  used  for  PRINT-filename:  otherwise,  prompt 
for  filename  OR  use  command-line  argument. 

DEFAULT:  FPRINT  =  blank 
UNITS:  none 


HMAX  Maxim.um  physical  altitude  f/max 

Maximum  physical  region  of  SSFPE  vertical  grid:  the  transmitter  altitude.  c0.  and 
all  output  receiver  altitudes.  cr.  must  be  less  than  Ffmax.  The  input  refractivity 
profiles  N( r)  must  all  extend  to  at  least  this  altitude.  In  addition.  #max  should 
exceed  the  height  of  the  highest  elevated  refractivity  profile  duct  (if  any)  along  the 
propagation  track. 

HMAX  =  0.  compute  ffniax  based  upon  the  highest  input  refractivity  profile  point 
along  propagation  track 

HMAX  >  0.  use  as  specified 
DEFAULT:  HMAX  =  0 
UNITS:  m 
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HMIN  Minimum  surface,  elevation  Hmjn 

Minimum  altitude  in  SSFPE  grid  for  variable  terrain  case:  Hmm  must  be  <  mini¬ 
mum  terrain  elevation  point  along  propagation  track.  All  input  profiles  must  start 
at  or  below  this  altitude. 

DEFAULT:  HMIN  =  0 

UNITS:  meters  relative  to  mean  sea  level 


ICURVE  Earth  curvature  flag 

Earth-flattening  transformation  flag:  controls  whether  the  earth-flattening  trans¬ 
form  is  to  be  used.  If  specified,  then  the  refractive  index  profiles  and  transmit¬ 
ter/receiver  heights  are  modified  to  account  for  an  earth-flattened  coordinate  sys¬ 
tem. 

ICURVE  =  0.  flags  a  spherical  earth  run:  modify  the  refract ivity  profile.  .Y(  r ). 
and  transmit  ter/receiver  altitudes  to  account  for  earth  curvature.  If  modified 
refractivity  M(z.  r )  is  input,  then  ONLY  correct  the  transmitter/receiver  altitudes. 

ICURVE  =  1  flat-earth  run:  no  earth-curvature  corrections  are  applied  to  either  the 
refractivity  profile  or  the  transmit  ter /receiver  altitudes 

RANGE:  ICURVE  =  {0.1} 

DEFAULT:  ICURVE  =  0 
UNITS:  none 

IDEBUG  Debug  print  flag 

Controls  quantity  and  type  of  debug  output  from  run:  output  is  to  file  specified 
by  FPRINT  variable  or  to  second  filename  on  command  line. 

IDEBUG  =  —  1.  no  PRINT-file  output 
IDEBUG  =  0.  no  diagnostic  debug 

IDEBUG  =  1.  minimal  debug  output  -t-  tabular  dB  loss  table 
IDEBUG  >  1.  detailed  debug  output 
DEFAULT:  IDEBUG  =  1 
UNITS:  none 


IGRAPH  Screen  graphics  flag 

Controls  the  run-time  screen  graphics  output.  Plot  layout  is  determined  by  vari¬ 
ables  in  the  PEGRAF-type  namelist. 

IGRAPH  =  0.  no  run-time  screen  graphics 
IGRAPH  =  1.  run-time  screen  graphics 
DEFAULT:  IGRAPH  =  0 
UNITS:  none 
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ILOSS  PE  output  flag 

Controls  type  of  VTRPE  model  output  from  the  run.  Output  is  dB-loss  values  on 
a  range- altitude  grid  determined  by  the  ZR-array  and  by  the  output  ranges. 
ILOSS  =  1.  output  path  loss  PL  =  20  log(2A‘0J?/|F|) 

ILOSS  =  2.  output  pattern  propagation  factor  PF  =  20  log  \F\ 

DEFAULT:  ILOSS  =  1 
UMTS:  none 

IPLOT  PLOT- file  flag 

Controls  format  of  the  output  PLOT-file  specified  by  the  FPLOT  variable  or  by 
third  filename  on  the  command  line.  The  PLOT-file  contains  dB  values  on  range- 
altitude  grid. 

IPLOT  =  0.  no  PLOT-file  produced 
IPLOT  =  1.  ASCII  file:  packed  integer  format 
IPLOT  =  2.  EREPS  compatible  format 
IPLOT  =  3.  binary  output  file 
DEFAULT:  IPLOT  =  0 
UNITS:  none 

IPOLAR  Transmitter  polarization  flag 

Flag  specifying  the  transmitter  field  polarization. 

IPOLAR  =  0.  horizontal  polarization 
IPOLAR  =  1,  vertical  polarization 
DEFAULT:  IPOLAR  =  0 
UNITS:  none 

IRUNIT  Range  units  flag 

Flag  specifying  input  range  units;  will  apply  to  all  refractivity  profile  and  surface 
data  unless  reset. 

IRUNIT  =  0,  ranges  in  kilometers 
IRUNIT  =  1.  ranges  in  nautical  miles 
IRUNIT  =  2.  ranges  in  kiloyards 
IRUNIT  =  3.  ranges  in  meters 
DEFAULT:  IRUNIT  =  0 
UNITS:  none 

ISOURCE  Transmitter  antenna  radiation  pattern  flag 

Specifies  type  of  transmitter  antenna  radiation  pattern  used  to  initialize  SSFPE 
field. 

ISOURCE  =  0.  Gaussian  antenna  pattern 
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I  SOURCE  =  1,  compound  antenna  pattern  (must  also  specify  the  CBEAM  parameter) 
ISOURCE  =  3.  Kaiser-Bessel  pattern 
DEFAULT:  ISOURCE  =  0 
UNITS:  none 

IVOLUME  Atmospheric  attenuation  flag 

Atmospheric  volumetric  absorption  loss  flag. 

IVOLUME  =  0.  no  volume  absorption  loss 

IVOLUME  =  1.  add  volume  loss  as  specified  in  ALPHA  array 

DEFAULT:  =  0 

UNITS:  none 

IVTPE  Variable  terrain  flag 

Flag  specifying  variable  surface  boundary  condition  run.  Lsed  to  specify  either  a 
variable  terrain  elevation  or  variable  surface  dielectric  properties  run. 

IVTPE  =  0.  earth's  surface  modeled  as  smooth  and  flat 

IVTPE  =  1.  surface  dielectric  properties  and/or  terrain  elevation  vary  with  range 
from  transmitter:  user  MUST  specify  GROUND-type  inputs. 

DEFAULT:  IVTPE=  0 

UNITS:  none 

IZUNIT  Vertical  units  flag 

Controls  all  input  vertical  units  for  transmitter/receiver  altitudes  and  refractivity 
profile  altitudes. 

IZUNIT  =  0,  all  input  altitudes  in  meters 
IZUNIT  =  1.  all  input  altitudes  in  feet 
DEFAULT:  IZUNIT  =  0 
UNITS:  none 

MEGAHZ  Transmitter  frequency 

Specifies  VTRPE  run  frequency  in  mega-hertz. 

DEFAULT:  Note:  mandatory  input 
UNITS:  megahertz 

NFFT  FFT  exponent  size  N 

Controls  the  FFT  size  used  in  run:  FFT  size  .V  =  2^^.  The  maximum  FFT  size 
allowed  is  dictated  by  the  amount  of  available  memory  on  the  computer. 

NFFT  =  0.  flags  a  variable  FFT  size  run:  program  selects  initial  FFT  transform 
size  based  on  input  values  of  ZMAX ,  DBPASS  and  wave  number  power  spectrum:  the 
FFT  size  is  then  dynamically  adjusted  as  the  run  progresses  to  ensure  no  aliasing. 
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NFFT  >  0,  flags  fixed  FFT  run:  program  runs  at  fixed,  specified  transform  size. 

This  option  also  disables  all  internal  error  checking.  CAUTION:  if  NFFT  is  too 

small,  then  PE  output  may  be  incorrect  due  to  transform  aliasing  problems. 

RANGE:  7  <  NFFT  <  12  PC  version 

RANGE:  7  <  NFFT  <  10  UNIX  version 

DEFAULT:  NFFT  =  0 

UNITS:  none 

NRANGE  Number  of  output  ranges  Nr 

Controls  the  number  of  range  points  Nr  at  which  pathloss  or  propagation  factor 
values  will  be  written  to  output  files.  Used  in  lieu  of  RLAST  to  specify  the  number 
of  output  ranges. 

NRANGE  =  0.  model  computes  number  of  output  ranges  Nr  as 

Nr  =  max[l.  (i?niax  - 

user  must  specify  RFIRST,  RLAST  and  RNGINC  in  this  case 

NRANGE  >  0.  specified  number  of  output  ranges,  with  last  range  being 

^max  =  -^0  +  (*V  ~  1  )^/? 


DEFAULT:  NRANGE  =  0 
UNITS:  none 

PESTEP  PE  range  step  size  Ar 

PE  algorithm  range  step  size.  WARNING:  Do  not  specify  this  parameter  unless 
sure  of  results. 

PESTEP  =  0,  flags  a  variable  step  size  run:  PE  range  step  Ar  is  computed  based 
upon  the  local  split-step  truncation  error  estimates.  The  range  step  size  will  be 
dynamically  adjusted  within  DRMIN  <  Ax  <  DRMAX  to  minimize  the  total  error  in 
the  PE  calculation 

PESTEP  >  0.  run  at  fixed  specified  range  step.  Warning  forcing  the  range  step 
size  disables  all  error  checks 
DEFAULT:  PESTEP  =  0 
UNITS:  meters 

RFIRST  First  output  range  /?0 

Minimum  range  of  PE  model  output:  if  /?0  =  0.  then  first  output  will  be  at  PE 
range  step  ri  that  is  nearest  to  A#  =  RNGINC. 

RFIRST  =  0.  first  output  range  i?0  =  RNGINC 
RFIRST  >  0,  first  output  range  /?0  =  RFIRST 
DEFAULT:  RFIRST  =  0 
UNITS:  units  as  specified  by  IRUNIT 
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RLAST  Last  output  range  i?max 

Used  in  lieu  of  Xr  =  NRANGE  to  specify  the  maximum  range  of  PE  output.  The 
number  of  output  range  points  Nr  is  then  computed  as 

A  /?  =  (7?max  -  R0)/Ar  +  1  . 

RLAST  =  0.  then  last  output  range  computed  by  model  as 

■ftmax  =  #0  +  (^r  ~  1)  *  '■ 

user  must  specify  NRANGE  parameter  in  this  case 
RLAST  >  0.  last  output  range  /?max  =  RLAST 
DEFAULT:  RLAST  =  0 
UNITS:  units  as  specified  by  IRUNIT 

RNGINC  Output,  range  increment  A# 

Range  increment  for  model  dB-loss  output.  L'sed  in  conjunction  with  /?y  =  RFIRST 
and  ,Vr  =  NRANGE  inputs  to  determine  the  ranges  at  which  loss  will  be  output  to 
the  PLOT  and  PRINT  files.  The  actual  output  range  {r( }  are  the  nearest  SSFPE 
range  step  to  the  user  specified  equi-spaced  range  values  7?,  =  i?0  +(/  — i)A„.  ;  = 

1-2 . -Vr. 

DEFAULT:  Note:  mandatory  input. 

RANGE:  RNGINC  >  0 

UNITS:  units  as  specified  by  IRUNIT 

TELEV  Trans  nutter  terrain  elevation  /?() 

Terrain  elevation  (relative  to  mean  sea  level)  at  the  transmitter: 

DEFAULT:  TELEV  =  0 

UNITS:  units  as  specified  by  IZUNIT 

THBMAX  Maximum  PE  vertical  aperture  0max 

Maximum  PE  vertical  angular  aperture  1 /2-width:  used  ONLY  for  compound  type 
aperture  starter. 

RANGE:  0  <  THBMAX  <  90 
DEFAULT:  none 
UNITS:  deg 

TITLE  Run  title 

Optional  run  title  used  to  label  output:  quote-delimited  character  string 
RANGE:  string  length  <  80  characters 
DEFAULT:  TITLE  =  N0SC  Radio  PE  Model 
UNITS:  none 
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VO  Sponge  amplitude  1’0 

SSFPE  complex  absorber  amplitude.  WARNING:  Do  not  specify  this  parameter 
unless  sure  of  results. 

VO  =  0.  absorber  parameters  computed  by  model 
VO  >  0.  specified  absorber  amplitude 
RANGE:  VO  >  0 
DEFAULT:  =  VO  =  0 
UNITS:  none 

WVO  Sponge  scale  length  tr0 

Vertical  length  scale  for  complex  “sponge”  absorber.  WARNING:  Do  not  specify 
this  parameter  unless  sure  of  results. 

WVO  =  0.  computed  based  upon  frequency  and  transmitter  properties 

WVO  >  0.  specified  absorber  scale  length 

RANGE:  WVO  >  0 

DEFAULT:  WVO  =  0 

UNITS:  meters 

ZMAX  Max  PE  grid  altitude  Zmax 

Maximum  e-value  in  PE  grid:  Zmax  =  ffmax  +  absorber  width.  Used  to  set  the 
vertical  wave  number  mesh  increment  Ap  via  A p  =  7r/Zniax.  To  avoid  aliasing  of 
Lloyd's  mirror  surface  interference  pattern,  requires  Znax  >  2zq.  where  is  the 
transmitter  altitude. 

ZMAX  =  0.  set  Zmax  based  upon  run  parameters 
ZMAX  >  0.  use  as  specified 
DEFAULT:  =  0 
UNITS:  meters 

ZR  Receiver  altitudes 

Receiver  (i.e.,  target)  altitudes  array  ZR():  dB-loss  values  (either  PL  or  PF)  are 
interpolated,  at  each  output  range,  from  the  internal  high-density  vertical  PE 
mesh  onto  a  vertical  grid  specified  by  the  contents  of  the  ZR-array.  The  entries  in 
ZR()  need  not  be  equi-spaced.  and  are  self-counting  —  the  first  non-positive  entry 
terminates  the  input  scan.  Two  options  are  available  to  specify  the  actual  receiver 
altitudes: 

(1)  input  a  list  of  V.  <  1000  monotonically  increasing  receiver  altitudes: 

ZR  ()  =Cj.C2 . ~\ 

(2)  specify  the  receiver  altitudes  by  inputting  via  ZR  up  to  five  sets  of  ordered 
triplets  having  the  form: 

ZR  =  Cj .  -  A  !  •  d]  [.  c2-  ^2  '  '  1 
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The  / tli  triplet,  {rj.  —  iV/.d/}.  will  then  be  expanded  into  Nj  equi-spaced  alti¬ 
tudes  starting  at  :  =  q  with  a  vertical  spacing  of  d/.  The  sum  of  all  the  Nfs 
must  be  <  1000. 

RANGE:  ZR(t)  >  0.  i  <  1000 

DEFAULT:  Note:  mandatory  in-put 

UNITS:  input  units  as  specified  by  IZUNIT  flag 

ZTRANS  Transmitter  altitude  ~0 
RANGE:  ZTRANS  >  0 
DEFAL’LT:  Note.:  mandatory  input 
UNITS:  input  units  as  specified  by  IZUNIT 

Table  1  summarizes  the  variables  in  the  RPE-type  namelist  and  their  defaults. 
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Table  1.  RPE  namelist  variables. 


Item 

Description 

Typea 

Units 

Default. 

ALPHA 

Atmospheric  attenuation  array 

R 

dB/km 

0. 

BMWIDTH 

Transmitter  vertical  beamwidth 

R 

deg 

none 

BMELEV 

Transmitter  vertical  beam  tilt 

R 

deg 

0. 

CBEAM 

Compound  aperture  parameter 

R 

0. 

DBPASS 

p-space  spectrum  pass-band 

R 

dB 

50. 

DRMAX 

Maximum  PE  range  step  size 

R 

m 

1000. 

DRMIN 

Minimum  PE  range  step  size 

R 

m 

75. 

EPSURF 

Surface  dielectric  constant 

C 

(0.,0.) 

FPLOT 

PLOT-filename 

A 

blank 

FPRINT 

PRINT- filename 

A 

blank 

HMAX 

Maximum  physical  altitude 

R 

m 

0. 

HMIN 

Minimum  surface  elevation 

R 

m 

0. 

I CURVE 

Earth  curvature  flag 

I 

0 

IDEBUG 

Debug  print  flag 

I 

1 

IGRAPH 

Screen  graphics  flag 

I 

0 

ILOSS 

PE  output  type  flag 

I 

1 

IPLOT 

PLOT-file  flag 

I 

0 

IPOLAR 

Transmitter  polarization  flag 

I 

0 

IRUNIT 

Range  units  flag 

1 

0 

ISOURCE 

Transmitter  antenna  type  flag 

I 

0 

IVOLUME 

Atmospheric  attenuation  flag 

I 

0 

IVTPE 

Variable  terrain  flag 

I 

0 

IZUNIT 

Vertical  units  flag 

I 

0 

MEGAHZ 

Transmitter  frequency 

R 

MHz 

none 

NFFT 

FFT  exponent  size 

I 

0 

NRANGE 

Number  of  output  ranges 

I 

0 

PESTEP 

PE  range  step  size 

R 

m 

0 

RFIRST 

First  output  range 

R 

IRUNIT'' 

0. 

RLAST 

Last  output  range 

R 

IRUNIT'1 

0. 

RNGINC 

Output  range  increment 

R 

IRUNIT'1 

none 

TELEV 

Terrain  elevation  at  transmitter 

R 

IZUNIT" 

0. 

THBMAX 

Maximum  PE  vertical  aperture 

R 

deg 

0. 

TITLE 

Run  Title 

A 

NOSC  ...  • 

VO 

Sponge  amplitude 

R 

0. 

WO 

Sponge  scale  length 

R 

IT) 

0. 

ZMAX 

Max  PE  grid  altitude 

R 

m 

0. 

ZR 

Receiver  altitudes  array 

R 

IZUNIT" 

none 

ZTRANS 

Transmitter  altitude 

R 

IZUNIT' 

none 

"FORTRAN  data  type:  A=character.  C=complex.  I=integer.  R=rea!. 
6 1 "nit s  set  by  IRUNIT. 

'  Units  set  by  IZUNIT. 
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§5.0  Atmospheric  Inputs 


The  SECTOR-type  namelist  variables  contain  the  atmospheric  refractivity  profile  data 
along  the  propagation  track.  For  a  range- independent  (i.e.,  single  refractivity  profile)  run. 
only  one  SECTOR-type  namelist  is  required:  for  a  range-dependent,  multi-profile  profile 
run.  a  separate  SECTOR-type  namelist  is  required  for  each  input  refractivity  profile.  Each 
refractivity  profile  input  defines  an  environmental  sector  identified  with  a  great  circle  range 
from  the  transmitter.  For  a  multi-profile  run.  each  SECTOR  namelist  must  be  ordered  by 
increasing  range  from  the  transmitter. 

5.1  Refractivity  Profiles 

The  VTRPE  model  requires  a  specification  of  the  atmospheric  radio  refractivity  field  as  a 
function  of  altitude  c  and  range  r.  The  refractivity  field  is  defined  by  one  or  more  vertical 
profiles  of  refractivity  vs.  altitude  at  various  ranges  from  the  transmitter.  Each  refractivity 
profile  is  input  as  ordered  pairs  {zj.  Pt  =  P(Cj)}.  where  zt  is  the  /til  profile  altitude  and  P- 
is  the  corresponding  refractivity  value.  The  refractivity  profile  P(z)  is  modeled  by  layers 
within  which  the  profile  is  assumed  to  be  piecewise  linear  with  altitude. 

The  refractivity  profile  may  be  specified  either  in  terms  of  refractivity  (X-units)  or  in 
terms  of  refractive  modulus  (M- units).  Each  refractivity  profile  that  is  input  must  extend 
from  the  minimum  physical  altitude  Hmin  (as  specified  by  the  HMIN  variable  in  the  RPE 
namelist),  to  the  maximum  physical  altitude  i/max  (as  specified  by  the  HMAX  variable  in 
the  RPE  namelist).  Linear  extrapolation  is  then  used  to  extend  the  profile  from  ifmax  to 
the  upper  boundary  of  the  PE  grid.  Zmax.  In  addition,  the  profiles  should  extend  to  a 
sufficiently  high  altitude  to  encompass  any  ducts  that  may  be  present.  The  last  layer  in 
the  input  profile  must  contain  a  positive  vertical  gradient.  All  profiles  should  be  specified 
to  an  altitude  beyond  which  simple  linear  extrapolation  is  appropriate. 

Optionally,  the  refractivity  profile  data  may  be  stored  in  separate  files  in  an  EREPS 
compatible  format,  and  the  filenames  specified  in  the  SECTOR  namelist. 

5.2  Profile  Range  Interpolation 

For  a  multiple  profile  run.  the  input  profiles  are  combined  to  form  a  continuous  refrac¬ 
tivity  field  in  range  and  altitude.  This  is  accomplished  by  means  of  a  bi-variate  surface 
fitting  algorithm  that  uses  the  input  pairs  {zj.Pj}  to  construct  a  smooth  surface. 

Range  interpolation  between  profiles  is  controlled  by  specifying  the  points  on  adjacent 
profiles  (vertices)  which  define  quadrilateral  surface  interpolation  regions.  If  no  surface 
interpolation  is  desired,  then  simply  connect  the  first  and  last  points  in  each  profile.  Note 
that  for  range-dependent  input.  ALL  profiles  must  be  extended  to  a  common  maximum 
altitude. 

5.3  SECTOR  Namelist  Variable  Description 

The  SECTOR  namelist  variables  are  defined  below: 
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CONNECT  Profile  vertex  connection  altitude.*  array 

Array  of  refractivity  profile  vertex  connection  altitudes  in  adjacent  profiles  that  are 
used  to  define  the  quadrilateral  regions  (tesselation)  used  for  profile  range/altitude 
interpolation.  Used  ONLY  for  range-dependent  data  inputs,  and  for  the  second 
and  successive  data  sectors.  The  CONNECT-array  data  are  input  as  ordered  pairs 
of  profile  vertices:  {c,(; ).  c2(/)}.  i  =  1.2.  •••  ,  NPAIR.  where  ~,  (/)  is  the  altitude 
of  the  /th  vertex  connection  in  the  previous  profile,  and  c2(0  is  the  altitude  of 
the  /th  vertex  connection  in  the  current  profile.  The  points  chosen  for  vertex 
connection  pairs  should  typically  be  local  extrema  (i.e..  local  maxima  or  minima) 
in  the  profile. 

DEFAULT:  none 

UNITS:  units  as  specified  by  IZUNIT 
IPTYPE  Profile  data  type  flag 

Specifies  type  of  atmospheric  refractivity  profile  data  that  are  input  via  the  PRO¬ 
FILE  array. 

IPTYPE  =  1.  contents  of  PROFILE  are  refractivity  ,Y(~)  in  N-units 
IPTYPE  =  2.  contents  of  PROFILE  are  modified  refractivity  M{ z )  in  M-units 
IPTYPE  =  3.  contents  of  PROFILE  are  normalized  modified  refractivity.  M(  z )  —  A/0. 
(M-units)  where  A/0  is  the  reference  modified  refractivity  as  specified  by  the  MZERO 
parameter 

DEFAULT:  IPTYPE  =  3 
UNITS:  none 

IRUNIT  Profile  range  units  flag 

Specifies  input  range  units  of  the  refractivity  profile  data:  overrides  value  of  IRUNIT 
set  in  RPE-type  namelist. 

IRUNIT  =  0.  profile  ranges  in  kilometers 
IRUNIT  =  1,  profile  ranges  in  nautical  miles 
IRUNIT  =  2,  profile  ranges  in  kiloyards 
IRUNIT  =  3,  profile  ranges  in  meters 
DEFAULT:  set  in  RPE  namelist 
UNITS:  none 

IZUNIT  Profile  altitude  units  flag 

Specifies  input  altitude  units  in  refractivity  profile:  overrides  IZUNIT  value  set  in 
RPE-type  namelist. 

IZUNIT  =  0.  input  heights  in  meters 
IZUNIT  =  1.  input  heights  in  feet 
DEFALLT:  set  in  RPE  namelist 
UNITS:  none 
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MZERO  Reference  modified  refractivity  Mq 

Reference  value  for  modified  refractivity:  MUST  be  input  when  IPTYPE=3.  Nor¬ 
mally  set  *o  the  surface  modified  refractivity  value. 

DEFAULT:  MZERO  =  0 

UNITS:  M-units 


NPAIR  Number  of  CONNECT  •point « 

Number  of  vertex  altitude  pairs  specified  in  the  CONNECT  array:  ONLY  used  for 
range-dependent  profile  inputs.  Should  not  be  specified  for  the  first  environmental 
sector! 

RANGE:  NPAIR  >  2 
DEFAULT:  none 
UNITS:  none 


NVP  Number  of  profi.lt  point*  X 

Number  of  points  in  the  input  profile  array  PROFILE:  MUST  be  specified  for 
multiple- profile  inputs  in  lieu  of  self-counting. 

NVP  =  0.  self-counting  input:  the  number  of  profile  input  points  determined  by 
the  model  based  upon  last  non-zero  entry  in  PROFILEQ.  The  self-counting  option 
may  only  used  for  single  profile  runs.  For  multiple  profile  runs,  each  profile  input 
must  specify  NVP. 

NVP  >  0.  number  of  data  pairs  input  via  PROFILE-array 
DEFAULT:  NVP  =  0 
L'NITS:  none 


PFILE  Profile  filename 

Optional  name  of  datafile  containing  input  refractivity  profile  data:  used  in  lieu 
of  inputting  the  profile  via  PROFILE  namelist  variable.  Quote-delimited  character 
string  containing  a  valid  filename  with  optional  path. 

PFILE  =  blank,  profile  input  via  PROFILE-array 

PFILE  =  filename",  read  profile  from  indicated  file  in  EREPS  format 
DEFAULT:  PFILE  =  blank 

PRANGE  Profile  range 

Profile  great  circle  range  relative  to  the  transmitter:  used  only  for  multiple-profile 
input.  Mandatory  input  for  a  multiple- profile  run:  PRANGE=0  for  1st  profile 

DEFAL’LT:  none 

UNITS:  as  specified  by  IRUNIT  variable 
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PRO  FID  Pro  jilt  ID 

Optional  profile  ID  character  string:  quote-delimited  character  string  (maximum 
SO  characters). 

DEFAULT:  PROFID  =  ‘  ‘ 

UNITS:  none 

PROFILE  Profile  array  P(z) 

Refract ivity  profile  data  array;  input  consists  of  ordered  pairs  {c-.P-}.l  <  i  < 
50.  where  c-  is  the  ith  profile  altitude  relative  to  mean  sea  level  and  Pt  is  the 
corresponding  refractivity  in  units  as  specified  by  the  IPTYPE  variable.  The  vertical 
units  are  specified  by  the  IZUNIT  variable.  The  first  profile  point  MUST  Ire  at  or 
below  the  minimum  terrain  elevation  c  =  Hmm  (Hmm  =  HMIN  is  set  in  the  RPE 
namelist)  along  the  propagation  track  for  a  variable-terrain  run.  or  at  the  surface 
(c  =  0)  for  a  smooth- surface  run.  The  last  point  in  the  profile  array  should  extend 
to  a  sufficiently  high  altitude  to  encompass  any  ducts  present  in  the  profile,  or  in 
any  event  be  greater  than  c  =  HMAX. 

DEFAULT:  none 

UNITS:  altitude  units  by  IZUNIT.  refractivity  units  by  IPTYPE 


Table  2  summarizes  the  variables  in  the  SECTOR-type  namelist. 

Table-2:  SECTOR  Namelist  Variables 


Item 

Description 

Type" 

Units 

Default 

CONNECT 

Profile  vertex  connection  array 

R 

IZUNIT6 

0. 

IPTYPE 

Profile  data  type  flag 

I 

3 

IRUNIT 

Profile  range  units  flag 

I 

RPE 

IZUNIT 

Profile  altitude  units  flag 

I 

RPE 

MZERO 

Reference  modified  refractivity 

R 

M-units 

0 

NPAIR 

Number  of  CONNECT  points 

I 

none 

NVP 

Number  of  profile  points 

I 

0 

PFILE 

Optional  profile  filename 

A 

blank 

PRANGE 

Profile  range 

R 

IRUNTT 

0. 

PROFID 

Profile  ID 

A 

blank 

PROFILE 

Refractivity  profile  array 

R 

IZUNIT6 

IPTYPE(/ 

none 

"  FORTRAN  data  type:  A=character.  C= complex.  I=integer.  R=real. 
6Units  set  by  IZUNIT. 
rVertical  units  set  by  IRUNIT. 

^Profile  units  set  by  IPTYPE. 


§G.O  Surface  Inputs 


Surface  environmental  inputs  are  specified  in  the  GROUND  namelist  input.  It  contains 
information  specifying  the  terrain  elevation  and/or  the  surface  dielectric  properties  along 
the  propagation  track.  This  information  is  used  by  the  VTRPE  model  if  a  variable  terrain 
run  is  specified  by  means  of  the  IVTPE  parameter  in  the  RPE  namelist.  If  the  earth's 
surface  is  not  flat  or  if  finite  conductivity  surface  boundary  conditions  are  desired,  then 
GROUND-type  input  data  must  be  supplied.  This  includes  the  terrain  elevation  vs.  range 
profile  and  specification  of  the  complex  surface  dielectric  constant. 

6.1  Terrain  Elevation  Data 


For  propagation  over  variable  terrain,  the  user  must  specify  the  terrain  elevation  vs. 
range  profile  h(r).  The  surface-elevation  profile  may  be  specified  in  one  of  two  ways: 

•  the  terrain  elevation  profile  h(r)  is  input  as  ordered  pairs  {r-.  /> • )}  where  r;  is  the 
f th  range  and  ht  =  h(r;)  is  the  corresponding  terrain  elevation  relative  to  mean  sea 
level.  Piecewise  linear  interpolation  is  used  to  construct  a  continuous  (with  range) 
profile  /((r)- 

•  by  specification  of  the  parameters  of  an  analytical  “bump"  surface-elevation  model. 
The  specific  analytic  form  used  is  a  Lorentzian  function  with  specified  peak  and 
FWHM  (full- width  at  half  maximum). 

6.2  Surface  Dielectric  Data 


The  complex  surface  dielectric  constant  is  defined  to  be  piecewise  constant  with  range 
and  may  be  specified  in  two  general  ways. 

First,  if  the  actual  complex  dielectric  constant  is  known,  then  it  is  input  directly  into 
the  VTRPE  model.  This  is  accomplished  by  specifying  the  surface  permittivity  e j  (i.e..  the 
real  part  of  the  complex  dielectric  constant)  and  either  the  imaginary  part  of  the  dielectric 
constant.  £2.  or  the  loss  tangent.  <*>.  where  tail 6  =  c2 /?\-  This  is  the  prefered  approach, 
provided  the  user  has  sufficient  knowledge  about  the  dielectric  properties  of  the  surface. 

Alternatively,  if  knowledge  of  the  complex  surface  dielectric  is  not  available,  then  the 
VTRPE  model  will  compute  surface  dielectric  constant  data  internally  using  semiempir- 
ical  models.  These  models  of  the  dielectric  data  are  based  upon  environmental  surface 
properties  such  as  temperature  and  salinity  that  are  commonly  available.  The  particular 
semiempirical  model  used  depends  upon  whether  the  surface  is  land  or  ocean. 

6.2.1  Sea  Water  Model 


If  the  surface  is  the  ocean,  then  a  semiei  ical  mode1  is  used  that  is  parameterized  by 
frequency,  sea  surface  temperature  and  sea  surface  salinity. 

For  a  good  conductor  like  sea  water,  the  frequency-dependent  dielectric  constant  is 
modeled  by  the  Debye  expression21 
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where  str  =  4.9  is  the  far- infrared  value  of  the  pure- water  dielectric  constant.  is  the 
static  dielectric  constant  of  sea  water  (i.e..  the  zero  frequency  limit),  r  is  the  relaxation 
time,  and  a  is  the  ionic  conductivity.  For  sea  water.  t0.  a.  and  r  are  all  functions  of  the 
teineperature  (T)  and  salinity  (s)  of  the  water. 

Klein  and  Swift22  have  analysed  L-band  experimental  data  to  derive  the  following  em¬ 
pirical  formulas  which  are  used  internally  by  the  VTRPE  model: 

c0  =  ^S7.134  —  0  1949 T  —  0.01276T2  +  2.491  x  10"4T3)  a(T.s). 

where 

a(T.s)  =  1.0  +  1.613  x  10_5sJ  -  3.656  x  HT3*  +  3.21  x  10~5s2  -  4.232  x  10-V. 

t=  (l.76Sx  HT11  -  6.086  x  10_13r  +  1.104  x  10~14r2- 8.111  a’O-17!3)  b(T.s), 
where 

6(T.s)  =  1.0  +  2.282  x  HT+sT  -  7.638  x  1(T4*  -  7.760  x  10“V2  +  1.105  x  lO-8*3. 
and 

a  =  <t(25.s)6_JA,  A  =  25  -  T. 

with 

3  =  d(A.s)  =  2.033  x  10-2  +  1.266  x  10_4A  + 2.464  x  10_6A2 
-.s  1.849  x  10-5  -  2.551  x  10"rA(  1.0 -0.1A)  . 

and 

<7(25.*)  =  *  (0.182521  -  1.46192  x  10-3*  +  2.09324  x  10-5*2  -  1.28205  x  10~7*3)  . 

In  the  above  formulas.  T  is  the  water  temeperature  in  °C’.  and  ,s  is  the  salinity  in  psu. 

6.2.2  Soil  Model 

For  propagation  over  ground,  a  wet  soil  model  developed  by  Hallikainen  et  al  is  used.23 
This  model  treats  the  ground  as  a  mixture  of  soil  particles,  air  voids,  and  liquid  water. 
The  water  contained  in  the  soil  is  divided  into  two  fractions:  ( 1)  bound  water,  and  (2)  free 
water.  Based  on  experimental  measurements  between  1.4  and  18  GHz.  the  dielectric  data 
at  each  frequency  are  fit  by  polynomial  regression  to  an  expression  of  the  form 

t  =  o0  +  AjS  +  a2C  +  (b0  +  b]S  +  b2C)mv  +  (c0  +  Cj  S  +  c2C)m2  , 

where  5  is  the  percentage  of  sand  in  the  soil.  C  is  the  percentage  of  clay  and  m  is  the 
soil  volumetric  moisture  content. 


6.3  GROUND  Input  Data  Stream 


The  user  may  input  surface  terrain  and  dielectric  data  in  one  of  six  ways: 

(1)  by  direct  input  of  the  complex  surface  dielectric  constant  es  —  € j  +  iz-2- 

(2)  by  input  of  the  sea  surface  temperature  and  salinity  (marine  only):  an  internal  com¬ 
putation  of  the  complex  sea  water  dielectric  constant  e  is  done  using  semiempirical 
models: 

(3)  by  input  of  the  surface  soil  type  and  moisture  content:  an  internal  computation  of 
the  complex  soil  dielectric  constant  via  semiempirical  formulas: 

(4)  by  input  of  the  surface  permittivity  5|  and  conductivity  a: 

(5)  by  input  of  the  surface  permittivity  and  loss  tangent  tan <5  =  or 

(6)  by  the  analytical  bump  model. 

The  actual  format  of  GROUND-type  inputs  is  a  pseudo  namelist  format.  The  input  file 
run  stream  is  scanned  for  the  starting  control  string.  &GR0UND.  and  the  terminating  control 
string.  "&END  ".  Each  control  string  must  be  on  a  separate  input  record,  start  in  column 
two  or  later,  and  have  no  embedded  blanks.  The  actual  surface  data  are  found  between 
these  two  control  strings  as  a  series  of  individual  data  records,  one  for  each  range,  preceded 
by  a  header  record.  The  input  format  of  the  header  record  and  data  records  is  standard 
FORTRAN  list-directed  input.  If  default  inputs  are  assumed,  then  commas  should  be  used 
to  specify  the  default  inputs.  All  header  and  data  records  should  be  terminated  with  a 
slash  (/).  The  GROUND  input  data  stream  is  organized  as  follows: 

&GR0UND  4=  starting  control  strn.g 

header  record 
data  record  1 

data  record  2  <=  input  data  records 


&END  <=  terminating  control  string 

6.3.1  Header  Record  Format 

The  header  input  data  record  has  the  format: 

IGTYPE,  IZUNIT,  IRUNIT 

The  first  item  in  the  header  record  is  the  data  type  flag  IGTYPE.  This  determines  the  type 
of  input  data  records  which  follow.  Each  input  data  record  must  have  the  same  format. 
Depending  upon  the  specification  of  the  IGTYPE  flag,  the  following  types  of  input  data 
records  are  allowed: 

(1)  IGTYPE=0:  range,  surface  elevation,  and  complex  surface  dielectric  constant  are 
input : 

(2)  IGTYPE=1:  range,  surface  elevation,  surface  permittivity,  and  surface  conductivity 
are  input: 
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(3)  IGTYPE=2:  range,  surface  elevation,  surface  permittivity,  and  the  surface  loss- 
tangent  are  input; 

(4)  IGTYPE=3:  range,  sea  surface  temperature,  and  salinity  are  input  (water  only) 

(5)  IGTYPE=4:  peak  range,  peak  height,  peak  width,  and  surface  dielectric  constant 
are  input  (bump  model); 

(6)  IGTYPE=5:  range,  surface  elevation,  and  the  soil  clay/sand/moisture  content  are 
input  ( semiempirical  soil  dielectric  model); 

IRUNIT  =  input  range  units  flag  (same  as  RPE  namelist) 

IZUNIT  =  input  z-units  flag  (same  as  RPE  namelist) 

6.3.2  Data  Record  Format 

Each  input  data  record  consists  of  a  range  value  ri  followed  by  a  list  of  parameters 
{ paranij }  that  specify  the  surface  properties  and  a  trailing  slash: 

range,  param j  [.param.2,parm 3....]  / 

The  format  is  FORTRAN  list-directed  input,  terminated  by  a  slash  (/). 

The  actual  order  and  contents  of  the  items  in  the  data  record  depend  upon  the  IGTYPE 
flag  in  the  header  record.  The  first  item  in  the  data  record  is  always  the  range.  Follow¬ 
ing  the  range  item  are  one  or  more  items  in  the  order  specified  in  the  following  list  as 
determined  by  the  IGTYPE  flag: 

IGTYPE  =  0  ri,hi.sl.s2: 

IGTYPE  =  1  r,-.  hh  cj.  a: 

IGTYPE  =  2  r,-, /i,-.  ij,  tan  6; 

IGTYPE  =  3  rpT.s; 

IGTYPE  =  4  Tj.  peak-height,  peak-width.  tj,t2  : 

IGTYPE  =  5  r{.  h{.  C.  S.  mv. 

The  input  variables  have  the  following  meaning: 

C  =  soil  clay  textural  component  by  weight  ( (A  ) 
a  —  surface  conductivity  (Siemens/meter) 
hi  =  surface  elevation  relative  to  mean  sea  level 
e2  =  imaginary  part  of  dielectric  constant 
tan  6  =  dielectric  loss  tangent  =  z2le\ 

mv  =  soil  volumetric  moisture  content  ^cm3/cm3j 
c  j  =  relative  dielectric  permittivity 
s  =  surface  salinity  (psu) 

5  =  soil  sand  textural  component  by  weight  (% ) 

T  =  surface  temperature  (centigrade) 


39 


§7.0  Output  Files 


There  are  three  possible  output  files  from  the  YTRPE  model: 

(1)  An  ASCII  forinatted-type  file  containing  a  diagnostic  listing  controlled  by  the  RPE 
namelist  input  variable  IDEBUG. 

(2)  An  optional  ASC’II-type  plot  file  controlled  by  the  RPE  namelist  input  variable 
IPLOT.  This  plot  output  file  is  intended  to  be  input  to  a  suitable  plotting  pro¬ 
gram  and  contains  tabular  values  of  path  loss  ( PL )  or  propagation  factor  ( PF) 
on  a  range/altitude  matrix  as  specified  by  inputs  in  the  RPE  namelist.  (See  the 
description  at  the  end  of  this  section.) 

(3)  An  optional  BIX  ARY- type  file  containing  path  loss  or  propagation  factor  on  a 
range /altitude  matrix. 

7.1  Plot  Output  File:  ASCII  Format 

The  plot  output  file  is  an  ASCII-type  file  consisting  of  a  header  block  followed  by  colum¬ 
nar  values  of  range  and  dB  loss  values  (i.e..  PL  or  PF).  The  header  block  is  comprised 
of  the  following  records  formed  by  five  FORTRAN  format ted-type  records  which  may  be 
read  using  list  directed  INPUT/OUTPUT  as  and  is  shown  in  figure  2. 
logical  record  1  run_title 
logical  record  2  program  version 
logical  record  3  input  Jilespec 
logical  record  4  integer  run  parameters 
logical  record  5  real  run  parameters 
logical  record  6  refractivity  profile  data 
logical  record  7  receiver  altitudes 

Figure  2.  Output  file  header. 

The  header  block  is  comprised  of  the  following  records  formed  by  five  FORTRAN 
formatted-type  records  which  may  be  read  using  list  directed  INPUT/OUTPUT  as: 

record  1:  run-title : 

quote  delimited  FORTRAN  character  string  containing  the  title  of  the  run 
(<  80  characters). 

record  2:  version-id  date -time -group  datajormat 

three  quote  delimited  FORTRAN  character  strings  in  the  order  indicated 
where 

•  version-id  is  the  VTRPE  program  version  identifica  .jn. 

•  date -time -group  is  the  date  and  time  of  the  run. 

•  data-f  ormat  is  the  FORTRAN  format  of  subsequent  data, 
record  3  input -file  spec 

full  name  of  file  containing  run  inputs, 
record  4  A  item  j  itf  n?2  •  •  •  item  y 

integer  parameters. when'  A  is  the  total  number  of  items  which  follow  (cur¬ 
rently.  -V  =  14)  and 
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item j  =  input  range  units  flag  IRUNIT 
item-2  =  input  altitude  units  flag  IZUNIT 
item 3  =  number  of  receiver  altitudes  AT 
item 4  =  number  of  output  ranges  Xr 
item 5  =  number  of  points  in  first  refract ivity  profile.  Xp 
item 6  =  profile  data  type  flag  IPTYPE 
item-  =  transmitter  antenna  type  flag  ISOURCE 
itemg  =  transmitter  polarization  flag  IPOLAR 
item$  =  type  of  absorber 
/fern jo  =  type  of  dB-loss  output  flag  ILOSS 
item jj  =  FFT  exponent  size  used  NFFT 
item  12  =  scale  factor  for  loss  values  ISCALE 
item} 3  =  number  of  ""triplets”  in  receiver  altitudes  array 
item 14  =  number  of  items  in  floating  point  record.  Xf[oa( 
record  5  item  j  item-2  ' ' "  item  \ 

floating  point  parameter  values,  where 
item j  =  transmitter  altitude  (m) 
it  cm2  =  transmitter  frequency  (Hz) 
item^  =  first  output  range  (km) 
item4  =  last  output  range  (km) 
item-g  —  output  range  increment  (km) 
itemg  —  transmitter  vertical  mainlobe  tilt  (deg) 
item-  =  transmitter  vertical  beamwidth  (deg) 
itemg  =  initial  PE  range  step  (m) 
item 9  =  maximum  altitude  in  PE  grid  (m).  Zmar 
itemiQ  =  complex  absorber  amplitude,  V0 
item}}  =  complex  absorber  scale  length. 
item j2  =  FFT  p-space  pass-band  limit.  DBPASS 

Following  the  floating  point  parameters  record,  are  one  or  more  records  which  contain 
Ap  profile  data  pairs  corresponding  to  the  first  input  refractivitv  profile.  Each  profile 
data  record  contains  up  to  five  pairs  of  altitude.  and  corresponding  refractivitv.  P,. 
{cj  P}  £2^2  '  }'  stafting  at  the  surface.  The  type  of  refractivitv  data  is  indicated  by 
itemg  of  the  integer  parameters  record. 

Following  the  profile  data  logical  record,  is  the  receiver  altitudes  logical  record  consisting 
of  one  or  more  physical  records  containing  the  AT  receiver  altitudes.  This  concludes  the 
header  block. 

After  the  header  block,  there  are  A'P  logical  range  records,  each  of  which  contain  A  Z 
scaled  (by  scale  factor  ISCALE)  integer  dB-loss  values.  Each  logical  range  record  comprises 
one  or  more  physical  records  and  has  the  format: 

range  item}  iteni2  •••  item 


41 


rang (  =  floating  point  range1  in  units  specified  by  IRUNIT 

ittnij  =  integer  scaled  dB-loss  value  corresponding  to  ith  receiver  altitude 

The  actual  FORTRAN  format  used  to  write  the  range  data  record  is  contained  in  the 
CHARACTER  string  data  ^formal  contained  in  logical  record  2  of  the  header. 

7.2  Plot  Output  File:  EREPS  Format 

If  IPL0T=2  is  specified  in  the  RPE  namelist,  then  the  output  plot  file  will  be  in  a  format 
that  may  be  input  into  the  EREPS  model.  Consult  the  EREPS  documentation  for  the 
exact  format. 

7.3  Plot  Output  File:  Binary  Format 

If  IPL0T=3  is  specified  in  the  RPE  namelist,  then  the  output  plot  file  will  be  a  FORTRAN 
binary-type  file.  The  file  will  consist  of  ,Vr  logical  records,  each  having  the  form 

range  j  (IB  j  dB-^  ■■■  dB  y 

where  range  (  is  the  PE  output  range,  and  the  {dBj}  are  the  model  output  db-loss  values. 


§S.()  Examples 


The  following  examples  of  VTRPK  input  run  streams  illustrate  various  options  previously 
described. 

8.1  Test  Case  1:  Evaporation  Duct 

Test  case  1  is  a  range-independent.  28-meter  evaporation- duct  profile,  with  smooth 
surface  boundary  conditions.  The  title  of  the  run  is  TEST  case  1:  28m  Evap  Duct.  A 

400  MHz  (MEGAHZ=400)  transmitter  is  located  at  a  height  of  25  meters  (ZTRANS=25)  above 
the  ocean  surface.  The  transmitter  has  a  Gaussian-type  beam  pattern  ( ISQURCE=0)  with  a 
vertical  beamwidth  of  2.5  degrees  ( BMWIDTH=2 . 5 )  and  is  steered  horizontally  (BMELEV=0). 
There  are  10  receivers  ( i.e. .  targets)  equispaced  at  0.5-meter  intervals  starting  at  1  meter 
(ZR=1.0,  -10.  ,  .5).  The  maximum  physical  altitude  in  the  calculation  is  309.5  meters 
(HMAX=309 .5).  Output  will  consist  of  dB  values  of  path  loss  (IL0SS=1 )  starting  at  the  origin 
(RFIRST*0).  spaced  at  0.25  km  (RNGINC=.25)  increments  (or  the  closest  actual  PE  step) 
out  to  a  maximum  range  of  50  km  (RLAST=50).  No  plot-fill-  will  be  produced  (IPL0T=0). 
The  refractivity  profile  is  input  in  scaled  modified  refract ivity  format  (IPTYPE=3).  with  a 
reference  modified  refractivity  value  of  309.5  M-units  (MZER0=309 . 5 ).  The  profile  consists 
of  23  height /refractivity  pairs,  with  the  vertical  units  being  meters. 

The  following  lines  show  the  input  file  corresponding  to  this  test  case.  Note  that  the 
namelist  identifier  strings  &RPE  and  &SECT0R  MUST  start  in  column  2. 

&RPE 

TITLE=  ’TEST  case  1:  28m  Evap  Duct’, 

IPL0T*  0, 

IS0URCE=  0, 

IL0SS=  1, 

IDEBUG*1 , 

MEGAHZ=  400.0, 

ZTRANS=  25.0, 

ZR=  1.0,  -10.0,  .5, 

BMWIDTH=  2.5,  BMELEV=  0.0, 

RFIRST=  0.,  RNGINC=  0.25,  RLAST=  50.0, 

HMAX  =  309.5, 
ftEND 

ASECT0R  IPTYPE*  3,  MZER0*  339.0, 

PROFILE* 

0.0,  0.0,  0.2,  -12.83,  0.316, 

0.794,  -17.08,  1.259,  -18.46,  1.995, 

5.012,  -22.31,  7.943,  -23.38,  12.589, 


-14.24, 

-19.81, 

-24.24, 


0.501, 

3.162, 

19.953, 


-15.67, 

-21.1, 

-24.76, 


43 


24.84,  39.811,  -24.48,  50.119,  -23.91,  63.096,  -23.01, 
21.69,  100.0,  -19.84,  125.893,  -17.32,  158.489,  -13.97, 
-9.56,  209.526,  -8.38,  309.526,  3.42, 


8.2  Test  Case  2:  Guadalupe  Island 

Test  case  2  is  an  example  of  a  range-dependent  run  over  water.  The  environmental 
data  correspond  to  the  Gaudalupe  Island  experiment  conducted  off  the  coast  of  South¬ 
ern  California  in  194S.  The  experiment  consisted  of  a  beacon-equipped  airplane  flying  a 
"saw  tooth  profile  away  from  receivers  mounted  on  a  bluff  overlooking  the  ocean.  Using 
reciprocity,  the  received  field  patterns  may  be  computed  by  reversing  the  position  of  the 
receiver  and  source. 

The  transmitter  properties  are  thus: 

( 1 )  frequency  of  200  MHz: 

(2)  Gaussian-type  antenna  pattern  having  a  2  degree  vertical  beamwidth  and  steered 
up  2  degrees: 

(3)  elevation  of  150  meters. 

The  range  units  are  nautical  miles,  with  path  loss  PL  output  every  2  nmi.  Ten  equi- 
spaced  receiver  altitudes  are  specified,  spaced  20  meters  apart  starting  at  10  meters.  Five 
refract ivity  profiles  are  input,  spaced  from  0  to  200  nmi  from  the  transmitter.  The  profile 
altitude  units  are  in  feet,  with  refract  ivity  in  N  units  being  specified.  The  first  profile  has 
28  points,  starting  at  the  sea  surface  and  extending  to  3500  feet.  Range  interpolation  is 
defined  between  the  first  profile  at  range  =  Omni  and  the  second  profile  at  range  =  SO  nmi 
by  the  specification  of  7  vertex  connection  pairs  that  determine  the  quadrilateral  bivariate 
surface  interpolation  regions.  The  tessellation  is  defined  as: 

( 1 )  between  r  =  0  on  profile  1  and  :  =  0  on  profile  2 

(2)  between  r  =  500  on  profile  1  and  r  =  500  on  profile  2 

(3)  between  r  =  500  on  profile  1  and  r  =  700  on  profile  2 

(4)  between  c  =  800  on  profile  1  and  r  =  1050  on  profile  2 

(5)  between  r  =  1500  on  profile  1  and  c  =  1500  on  profile  2 

(0)  between  r  =  1G00  on  profile  1  and  c  =  1600  on  profile  2 

(7)  between  z  =  3500  on  profile  1  and  r  =  3500  on  profile  2 

ftRPE 

TITLE=  '"Guadalupe  Island  Experiment’, 

MEGAHZ=  200., 

BMELEV=  2 . , 

BMWIDTH=  2. , 

ZR=  10. ,  -10. ,  20. , 

ZTRANS=  150., 


28.0, 
79.433,  - 
199.526, 
&END 
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ZMAX=1500  .  , 

IDEBUG=1 , 

IPL0T=2 , 

ILOSS=  1, 

NRANGE=100 , 

RNGINC3  2.  , 

IRUNIT  =  1, 

RFIRST=  0 .  , 

&END 

&SECTOR 

PROFID3’ Guadalupe  Is  Exp:  Modified  profile  #1’,  PRANGE=0. , 
IPTYPE=1 ,  IRUNIT=1 ,  IZUNIT=1 ,  NVP=  28, 

PROFILE3 


0. 

,  336 .  , 

200. 

,  335.  , 

300. 

,  334., 

400. 

,  334., 

500. ,  332. , 

550. 

,  329., 

600. 

,  326., 

650. 

,  314., 

700. 

,  302., 

800.,  285., 

850. 

,  284., 

900. 

,  282., 

1000 . 

,  279., 

1200. 

,  276., 

1300.,  276., 

1400. 

,  276., 

1500. 

,  275., 

1600. 

,  273., 

1800. 

,  272., 

2000.,  272., 

2200. 

,  271., 

2400. 

,  269.  , 

2600. 

,  267., 

2800. 

,  266., 

3000.  ,  266.  , 

3200. 

,  266., 

3400. 

,  264., 

3500. 

,  263., 

&END 

ASECTOR 

PROFID3’ Guadalupe  Is  Exp: 

profile  #2’ , 

PRANGE=80 . , 

NVP=28 , 

PROFILE3 

0.  , 

337.  , 

200.  , 

335.  , 

400.  , 

333.  , 

500.  , 

332.  , 

700.,  330., 

800.  , 

324.  , 

900.  , 

309.  , 

950.  , 

304.  , 

1000.  , 

296.  , 

1050. ,  293. , 

1100.  , 

291.  , 

1200. , 

290.  , 

1300. , 

287.  , 

1400. , 

286.  , 

1500. ,  282. , 

1550.  , 

275.  , 

1600. , 

273.  , 

1700. , 

273.  , 

1800. , 

272.  , 

2000. ,  272. , 

2200.  , 

271., 

2400. , 

269.  , 

2600. , 

267.  , 

2800. , 

266.  , 

3000. ,  266. , 

3200 .  , 

266.  , 

3400 . , 

264.  , 

3500. , 

263.  , 

NPAIR=7 , 

CONNECT3  0 . ,  0 .  , 

500. ,  500. , 

500. ,  700. , 

800.,  1050., 

1500.,  1500., 

1600.,  1600., 

3500. ,  3500. 

ftEND 

ASECTOR 


PROFID3 ’Guadalupe  Is  Exp: 

PROFILE3 

profile  #3’ , 

PRANGE=120 .  , 

NVP=28 

J 

0.  , 

334.  , 

200. ,  334. , 

300.  , 

333.  , 

400.  , 

333.  , 

600.  , 

331.  , 

700.  , 

330.  , 

800. ,  328. , 

1000.  , 

327.  , 

1100.  , 

327.  , 

1200.  , 

326.  , 

1250.  , 

319.  , 

1300. ,  312. , 

1400.  , 

295.  , 

1450.  , 

291.  , 

1500.  , 

286.  , 

1550.  , 

284.  , 

1600. ,  283. , 

1800.  , 

281 .  , 

1900.  , 

276.  , 

2000.  , 

271.  , 

2100.  , 

3000.  , 

270.  , 

265.  , 

2200. ,  269.  , 

3300. ,  263. , 

2400.  , 

3500.  , 

269.  , 

263.  , 

2600.  , 

269.  , 

2800. , 

266.  , 

NPAIR37 , 

CONNECT3  0.,  0., 

700. ,  700. , 

700.,  1200., 

1050.,  1500., 

1500.,  1800., 

2000.,  2000., 

3500. ,  3500. 

&END 

&SECTOR  PROFID3’ Guadalupe  Is  Exp:  profile  #4’, 


PRANGE=160 . 

,  NVP=24 , 

PROFILE3 

0 . ,  333 .  , 

50. ,  332.  , 

300. ,  330. , 

400. ,  329. , 

600. ,  327. , 

800. ,  325.  , 

1000.  ,  324.  , 

1150. ,  323. , 

1200. ,  320.  , 

1400. ,  305. , 

1500.  ,  297.  , 

1600. ,  292.  , 

1650. ,  289. , 

1700.,  287., 

1800. ,  282. , 

2000.  ,  273.  , 

2200.,  267., 

2400. ,  271. , 

2500. ,  273. , 

2600. ,  272. , 

2800.  ,  267.  , 

3000. ,  265.  , 

3300. ,  265. , 

3500. ,  261. , 

NPAIR=6 , 

CONNECT3  0.,  0.,  1000.,  1000.,  1200.,  1150.,  1500.,  1600., 

2000. ,  2000. ,  3500. ,  3500. 

&END 

ASECTOR  PROFID3’ Guadalupe  Is  Exp:  profile  #5’, 

PRANGE=200 . ,  NVP=24, 

PROFILE3 


0. ,  331. 

,  200. 

,  329. 

300. 

,  327., 

400. 

,  327., 

600.  , 

324.  , 

800. ,  321. 

,  1000. 

,  320. 

,  1200. 

,  318., 

1400. 

,  315., 

1600. , 

307.  , 

1700. ,  299. 

,  1750. 

,  301. 

,  1850. 

,  290., 

2000. 

,  280., 

2200. , 

274.  , 

2400  .  ,  270  . 

,  2500. 

,  273. 

,  2600. 

,  274., 

2700. 

,  278., 

2800. , 

265.  , 

2900.  ,  261. 

,  3100. 

,  264. 

,  3300 . 

,  263., 

3500. 

,  262., 

NPAIR=8 , 

CONNECT3 

0. 

0.  , 

1000. 

1000. , 

1200. 

1600. , 

1600. 

2000 
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2000.  2000.,  2200.  2200.,  2400.  2400.,  3500.  3500. 

&END 


8.3  Test  Case  3:  N.  Persian  Gulf 

This  rest  case  is  a  range  dependent,  variable  terrain  run  between  0.5  and  200  km.  The 
frequency  is  520  MHz.  with  the  transmitter  located  at  60  meters  and  1  receiver  at  10 
meters  altitude.  There  are  three  refract ivity  profiles  specified,  and  the  surface  properties 
are  specified  in  the  IGTYPE=0  format. 


&RPE 

TITLE=  *N.  Persian  Gulf  Test’, 
MEGAHZ=  520., 

NFFT=0 , 

BMELEV=  0 . ,  BMWIDTH=  1.6, 

ZR=  10. , 

ZTRANS=  60 . , 

DRMAX  =  500.  , 

DRMIN  =  75. , 

PESTEP=0 .  , 

HMAX=1000 .  , 

IDEBUG=2 , 

IPL0T=0 , 

IL0SS=  1, 

RFIRST= . 5 ,  RNGINC= . 5 ,  RLAST=200. 
IVTPE-1, 

&END 

&SECT0R 

PR0FID=’16m  Evap  Duct:  profile  #1’, 


PRANGE 

=0.  , 

IPTYPE=3 ,  MZERG: 

=339.  , 

NVP=  22, 

PROFILE* 

0.,  0., 

0.1,  -12.995, 

0.25, 

-14.839, 

0.50, 

-16.161, 

1.0, 

-17.485, 

2.0,  -18.746, 

4.0, 

-19.882, 

6.0, 

-20.443, 

8.0, 

-20.768, 

12.,  -21.08, 

14.0, 

-21.138, 

16.0, 

-21 . 155, 

18.0, 

-21.140, 

20.0,  -21.101, 

25.0, 

-20.923, 

50.  , 

-19.184, 

100.  , 

-14.320, 

125.0,  -11.64, 

150.0, 

-8.88, 

200.0, 

-3.206, 

1000.0, 

93.57, 

1500.,  155.26 

aEND 
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&SECTOR 

PROFID3’ Profile  #2’, 

PRANGE=100 . ,  IPTYPE=2 ,  MZER0=0 . ,  NVP=12, 

PROFILE= 


c. , 

376.0, 

50.  , 

363.0, 

100.  , 

339.0, 

125.  , 

325.0, 

150.  , 

315.0, 

200.  , 

320.  , 

250.  , 

321.  , 

300.  , 

319.  , 

350.  , 

318.  , 

380.  , 

325.  , 

500.  , 

339.  , 

1500.  , 

475.  , 

NPAIR=3 

> 

CONNECT 

= 

o., 

0., 

16.,  150 

1500.  , 

1500. 

&END 

&SECTOR 

PROFID= 

’Profile 

PRANGE3 

140.  , 

IPTYPE= 

2, 

MZER0=0 

•  > 

NVP3  9, 

PROFILE 

= 

0.  , 

275.0, 

150.  , 

315.0, 

200.  , 

320.  , 

250.  , 

321.  , 

300.  , 

319.  , 

350.  , 

318.  , 

380.  , 

325.  , 

500.  , 

339.  , 

1500.  , 

475.  , 

CONNECT3 
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0. , 

o.. 

150.  , 

150.  , 

1500.  , 

1500. 

&END 

&GR0UND 

0,  0, 

0 

o. , 

o. , 

70.  , 

50. 

/ 

98.  , 

o. , 

70.  , 

50. 

/ 

100.  , 

o. , 

70.  , 

50. 

/ 

101 .  , 

1.  , 

70.  , 

50. 

/ 

124.2 

1.5, 

70.  , 

50. 

/ 

126.0, 

2.0, 

24.  , 

32. 

/ 

129.8, 

60.  , 

10., 

1.6 

/ 

132.6, 

100.  , 

3., 

.02 

/ 

136.6 

200.  , 

3.  , 

.02 

/ 

138.2, 

160.  , 

3.  , 

.02 

/ 

142.6, 

100.  , 

3.  , 

.02 

/ 

157.4, 

100.  , 

3., 

.02 

/ 

158.6, 

200.  , 

3., 

.02 

/ 

159.4, 

300.  , 

3.  , 

o 

to 

/ 

160.4, 

400.  , 

3-, 

.02 

/ 

161.8, 

150.  , 

3., 

.02 

/ 

162.2, 

325.  , 

3., 

.02 

/ 

177.4, 

400.  , 

3.  , 

.02 

/ 

191.0, 

400.  , 

3.  , 

.02 

/ 

197.0, 

500.  , 

3., 

.02 

/ 

201 .  , 

600.  , 

3.  , 

.02 

/ 

202.  , 

500.  , 

3., 

.02 

/ 

212.6, 

600.  , 

3.  , 

.02 

/ 

250.6, 

600.  , 

3.  , 

.02 

/ 

ftEND 
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